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PREFACE 


This  report  is  the  second  of  a series  dealing  with  the  manipula- 
tion and  interpretation  of  Landsat  digital  data.  The  concepts  and  pro- 
cedures reported  herein  were  developed  under  the  In-House  Laboratory- 
Independent  Research  Program,  Project  i+Al6ll01A91D,  Task  02,  Work  Unit 
095,  "Feasibility  of  Using  Landsat  Spectral  Data  for  Acquisition  of 
Terrain  Information  for  Multiple  Purposes."  The  work  was  performed 
during  the  period  September  19T5-June  1976  by  personnel  of  the  Environ- 
mental Simulation  Branch  (ESB),  Environmental  Systems  Division  (ESD) , 
Mobility  and  Environmental  Systems  Laboratory  (MESL),  U.  S.  Army  Engi- 
neer Waterways  Experiment  Station  (WES),  under  the  direct  supervision 
of  Messrs.  H.  W.  West,  Project  Manager,  and  J.  K.  Stoll,  Chief,  ESB.  The 
study  was  under  the  general  supervision  of  Messrs.  B.  0.  Benn,  Chief, 

ESD,  and  W.  G.  Shockley,  Chief,  MESL.  Dr.  H.  Struve,  ESB,  was  responsi- 
ble for  the  development  of  the  interactive  procedure  for  classifying 
terrain  types  by  spectral  characteristics.  Special  acknowledgment  is 
made  to  Mr.  W.  E.  Grabau,  Special  Assistant,  MESL,  for  his  constructive 
criticisms  and  encouragement  throughout  the  study.  Dr.  Struve  and 
Messrs.  Grabau  and  West  prepared  the  report. 

COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE,  were  Directors  of 
the  WES  during  the  study  and  report  preparation.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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CONVERSION  FACTORS,  METRIC  (Sl)  TO  U.  S.  CUSTOMARY  AND 
U.  S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  he  converted  as  follows : 

Multiply To  Obtain 

Metric  (SI)  to  U.  S.  Customary 


micrometres 

3.937OOT  X 

H 

0 

inches 

millimetres 

O.O39370OT 

inches 

centimetres 

0.3937007 

inches 

metres 

5.399568  X 

10"^ 

miles  (U.  S. 

nautical ) 

metres 

6.213711  X 

10-"* 

miles  (U.  S. 

statute ) 

kilometres 

0.5399568 

miles  (U.  S. 

nautical ) 

kilometres 

0.6213711 

miles  (U.  S. 

statute ) 

square  kilometres 

21+7.105^ 

acres 

square  kilometres 
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square  miles 

(U.  S.  nautical) 

square  kilometres 

0.3861021 

square  miles 

(U.  S.  statute) 

hectares 

2.U71054 

acres 

U.  S.  Customary  to 

Metric 

(SI) 

feet 

0.3048 

metres 

degrees  ( angular ) 

0.01745329 

radians 
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ACQUISITION  OF  TERRAIN  INFORMATION  USING 
LANDSAT  MULTISPECTRAL  DATA 

AN  IHIERACTIVE  PROCEDURE  FOR  CLASSIFYING  TERRAIN 
TYPES  BY  SPECTRAL  CHARACTERISTICS 

PART  I ; INTRODUCTION 

Backp;round 

1.  Modern  inilitaxy  operations  require  vast  amounts  of  terrain  in- 
formation, both  for  long-range  strategic  planning  and  quick-response 
tactical  planning.  There  is  also  a host  of  more  prosaic  needs  for 
terrain  information,  such  as  long-range  planning  for  and  monitoring  of 
military  facilities  uses.  The  various  agencies  that  supply  the  U.  S. 

Ai’med  Forces  with  terrain  information  have  used  conventional  air-photo 

1 2 

interpretative  procedures  for  many  years.  ’ Unfortunately,  needs  have 
grown  more  rapidly  than  capabilities,  and  thus  a requirement  for  faster 
and  less  costly  ways  of  acquiring  the  necessary  tei’rain  information  is 
urgent.  One  possible  avenue  is  that  opened  by  the  Landsat  satellites 
and  their  proposed  follow-ons. 

2.  There  are  (1977)  two  satellites  (Landsat  1 and  Landsat  2)  in 
orbit  around  the  earth.  Each  satellite  circles  the  earth  every  103  min 
and  contains  a raultispectral  scanner  (MSS)  that  provides  radiance  mea- 
surements of  terrain  materials  for  four  spectral  bands  (O. 5-0.6,  0. 6-0.7, 
0.7-0. 8,  and  0. 8-1.1  pm*)  for  each  pixel  area  (“57  by  79  m)  of  tlie 
terrain  sm’face.  Spectral  data  for  an  area  on  the  ground  are  provided 
every  l8  days  by  each  orbiting  satellite,  and  the  two  satellites  are 
scheduled  in  ouch  a manner  that  actual  coverage  (i.e.  spectral  data)  is 


* A table  of  factors  for  converting  metric  (Sl)  luiits  of  measur'ement 
to  U.  S.  customary  units  and  U.  S.  customary  units  to  metric  (Gl) 
units  is  presented  on  page 
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obtained  every  9 days.*  Computer-compatible  tapes  (CCT's)  of  the  spec- 
tral data  and  imagery  made  from  the  tapes  are  made  available  by  the  Earth 
Resoui'ces  Observation  Systems  (EROS)  Data  Center**  approximately  6 weeks 
after  an  area  has  been  scanned  by  the  satellite  sensor  system. 

Objective  and  Scope 

3.  The  objective  of  the  research  described  herein  was  to  develop 
a rapid  and  semiautomated  method  of  classifying  Landsat  radiance  data 
in  terms  of  preselected  land-use  categories  or  terrain  types . The  pro- 
cedure is  intended  as  an  interim  solution  to  the  problem  of  mapping  very 

2 

large  areas  (in  excess  of  100,000  km  ) in  terms  of  relatively  crude 
categories  in  very  short  periods  of  time.  \ major  requirement  was  that 
the  procedure  operate  with  a minimal  amount  of  processing  of  primary 
Landsat  data.  Ideally,  it  would  operate  by  manipulating  Landsat  data  as 
it  is  recorded  by  the  satellite,  without  data  decompression,  calibra- 
tion, or  other  modification. 

Rationale  of  the  Procedure 


it . Any  semiautomated  procedure  must  be  somewhat  anagolous  to  the 
procedure(s)  used  by  human  interpreters;  therefore,  an  understanding  of 
the  various  methods  by  which  human  interpreters  identify  terrain  and 
land-use  information  from  images  is  necessary.  There  appear  to  be  three 
basic  techniques: 

a.  Identification  of  features  or  materials  by  spectral 

analysis.  This  process  assumes  prior  knowledge  of  the 
absorption  and  reflection  characteristics  of  the  featm'e 
for  which  a search  is  being  made.  Identification  is  made 


* At  the  time  this  report  was  being  published,  the  National  Aeronautics 
and  Space  Administration  (NASA)  was  rescheduling  the  Landsat  coverage 
intervals.  After  the  proposed  launch  in  1978  of  landsat  C,  Landsat  1 
will  follow  Landsat  2 six  days  later,  and  Landsat  C will  follow 
Landsat  2 nine  days  later. 

**  United  States  Department  of  the  Interior,  Geological  Survey,  EROS 
Data  Center,  Sioux  Falls,  S.  Dak.  57198. 
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by  matching  the  known  spectrum  with  the  spectrum  as 
revealed  by  the  sensor  system. 

b.  Identification  of  featm~es  or  materials  by  pattern 

recognition.  This  process  assumes  that  each  feature  or 
material  exhibits  a unique  pattern  of  colors  or  shades 
of  gray.  Recognition  is  achieved  when  a pattern  on  an 
image  can  be  correlated  with  an  object  or  material  on  the 
ground. 

c_.  Identification  of  features  or  materials  by  association. 
This  seems  to  be  the  major  process  used  by  the  human 
interpreter.  Despite  the  fact  that  it  is  widely  used,  it 
is  very  poorly  understood.  In  fact,  it  cannot  even  be 
defined  except  by  an  example:  A smooth  gray  area  around  a 
residential  building  is  almost  instantly  recognized  by  an 
interpreter  as  a lawn,  not  because  it  is  smooth  and  gray, 
but  because  it  is  positioned  in  a certain  way  with  respect 
to  the  building.  The  precise  same  shade  of  gray  associ- 
ated with  a building  identified  as  a factory  might  well  be 
interpreted  as  a parking  lot.  An  immense  amount  of  time 
is  devoted  by  educational  institutions  to  the  development 
of  the  photo  interpreter's  ability  to  make  associations, 
and  an  interpreter's  skill  is  largely  measured  by  the 
number  of  associations  he  can  reliably  make. 


5.  The  interactive  pi'ocedure  uses  the  developed  skill  of  a human 
interpreter;  the  human  makes  the  primary  interpretation  and  then  speci- 
fies to  a computer  program  those  combinations  of  spectral  values  that 
are  to  be  accepted  as  representing  the  terrain  types  of  interest  on  the 
ground.  The  basic  reason  for  exploiting  an  interactive  procedure  is 
that  the  human  interpreter  is  capable  of  m'>.king  subjective  judgments  on 
the  basis  of  all  three  of  the  techniques  described  above;  whereas  at  the 
present  level  of  development,  programmed  interpretation  can  be  based 
only  on  the  first  of  the  three  teclmiques  (i.e.  spectral  composition). 
The  advantage  is  tliat  classifications  can  sometimes  be  based  on  quite 
subtle  criteria,  such  as  those  tliat  ctinnot  be  revealed  or  defined  by 
rigorous  objective  analysis. 

6.  The  disadvantages,  however,  are  compelling.  Perhaps  tlie  most 
important  is  that  each  image  must  be  independently  subjected  to  the 
interpretative  process.  A classification  made  on  the  basis  of  an  image 
obtained  on  Monday  cannot  be  reliably  cai’ried  over  to  an  image  obtained 
on  Tuesday.  Nor  can  classifications  be  reliably  transposed  to  another 
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geographic  region,  even  if  close  by.  The  reasons  are  many.  The  most 
important  are  as  follows : 

a.  Spectral  characteristics  are  changed  by  differences  in 
atmospheric  composition. 

Spectral  characteristics  are  changed  by  differences  in 
reflectance-  geometry. 

c_.  Spectral  characteristics  are  changed  by  time-dependent 
variations  in  terrain  c Mditions. 

7.  The  skilled  human  interpreter  can  often,  somehow,  sort  out  and 
subjectively  compensate  for  all  of  these  things;  but  the  basic  digital 
data,  on  which  a computer  program  must  base  its  decisions,  may  be  quite 
dramatically  different  from  scene  to  scene.  Writing  a program  that  will 
automatically  make  the  necessary  adjustments  nd  accommodations  is  well 
beyond  current  capability.  However,  it  is  c..  f^ar  that  the  long-range 
goal  of  programmed  interpretation  is  a program  that  will  make  the  neces- 
sary adjustments.  A portion  of  such  a program  is  discussed  in  detail  in 
Reference  3. 
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PART  II;  DEVELOPMENT  OF  SEMIAUTOMATED  PROCEDURE 


Theory 

Primary  data  versus  processed  data 

8.  In  principle,  the  primary  record  stored  on  Landsat  CCT's  con- 
tains all  of  the  spectral  data  that  will  he  available.  The  primary  data 
may  he  modified  hy  calibrating  them  to  references  and  thus  converted  to 
true  radiance  values.  They  may  be  corrected  for  atmospheric  transmit- 
tance and  other  effects,  but  those  processes  merely  change  the  absolute 
values;  they  do  not  change  the  actual  informational  content.  Therefore, 
at  least  in  theory,  the  primary  OCT  record  is  as  useful  for  purposes  of 
terrain  classification  as  those  data  after  calibration  and  other  forms 
of  processing.  That  is,  if  "natural"  classes  exist,  they  will  be  in- 
herent in  the  primary  data  as  well  as  in  the  processed  data. 

9.  This  suggests  that  correlations  that  relate  natural  classes  of 
CCT  primary  data  with  conditions  on  the  ground  (i.e.  with  terrain  types) 
are  at  least  as  valid  as  those  obtained  with  processed  data.  If  this  is 
so,  then  no  advantages  are  to  be  gained  by  correcting  primary  data  in 
any  way  prior  to  classification.  Thus,  if  interpretation  is  to  be  based 
strictly  on  the  establishment  of  correlations  between  spectral  composi- 
tion of  pixels  (i.e.,  values  from  each  of  four  wavelength  bands,  if  the 
sensor  is  Landsat)  and  terrain  types  (not,  be  it  noted,  between  spectral 
composition  of  pixels  and  known  spectral  characteristics  of  the  terrain 
types),  then  the  primary  data  are  at  least  as  good  as  any  form  of  pro- 
cessed data. 

Determining  spectral  signatures 

10.  To  simplify  subsequent  discussion,  the  primary  Landsat  data 
defining  the  spectral  characteristics  of  a pixel  is  called  a "CCT  value 
set."  Thus,  a CCT  value  set  is  an  array  of  four  CCT  values,  one  for 
each  wavelength  band. 

11.  Let  it  be  assumed  that  the  CCT  value  set  reflected  from  a 
particular  terrain  type  always  consists  of  four  fixed  CCT  values  under 
any  given  set  of  conditions  (i.e.  of  atmospheric  transmittance  and 
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reflectance  geometry).  Now  suppose  that  a region  including  several  dif- 
ferent terrain  types,  each  always  exhibiting  a set  of  four  fixed  CCT 
values,  is  imaged  by  a multispectral  scanner,  such  as  that  aboard  Landsat 
If  the  CCT  values  of  any  one  of  the  four  spectral  bands  were  assembled 
into  a histogram  of  the  number  of  pixels  versus  the  CCT  values,  the 
histogram  would  consist  of  a set  of  discrete  bars,  each  representing  the 
number  of  pixels  exhibiting  a particular  CCT  value.  For  example,  a 
portion  of  a Landsat  scene  consisting  of  33,600  pixels  might  exhibit  the 
histogram  shown  in  Figure  1. 


TOTAL  SCENE  210  « 160  PIXELS 


Figure  1.  Idealized  histogram  of  one  wavelength  band 

12.  However,  experience  has  shown  that  CCT  value  sets  for  a given 
terraiii  type  tend  to  distribute  themselves  to  either  side  of  average 
values  rather  than  remain  fixed  values.  Therefore,  histogram  sets  of 
several  terrain  types  using  real  world  data  (Figure  2)  will  take  the 
form  of  a continuous  curve  rather  than  discrete  bars.  There  are  at 
least  three  major  sources  that  act  randomly  and  in  concert  to  cause 
the  CCT  values  to  fluctuate. 

a.  The  Landsat  scanner  has  an  internal  error  of  about  +2 
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togram  of  wavelength  bands  of  one  quarter  of  Landsat  scene 


percent  of  total  range.  Thus,  for  example,  a true 
radiance  value  of  30  may  be  recorded  as  any  value  between 
28.7  and  31.3.  Since  the  recorded  CCT  values  are  always 
integers,  the  values  become  29  and  31. 

b.  Local  variations  in  atmospheric  transmittance  and  terrain 
conditions  alter  the  amount  of  radiation  reaching  .the 
sensor. 


c . 


Pixels  that  fall  partly  in  one  terrain  type  and  partly  in 
another  will  yield  CCT  values  that  are  somewhere  between 
the  two.  Thus,  they  belong  to  no  "naturalistic"  classes. 


since  they  do  not  really  represent  a single  set  of  ground 
conditions.^  All  Landsat  images  contain  some  of  these 


"spurious"  pixels,  and  most  contain  very  many. 


13.  Despite  the  "smearing"  of  the  theoretical  values,  let  it  never- 
theless be  assumed  that  each  peak  represents  a CCT  value  characteristic 
of  a terrain  type.  It  might  also  be  assumed  that  a distinct  "shoulder," 
such  as  that  at  38  in  the  band-U  histogram  (Figure  2a),  would  represent 
a weakly  developed  peak,  possibly  because  the  area  covered  by  the  terrain 
type  is  small.  The  obvious  extension  of  this  notion  is  to  divide  the 
CCT  value  scale  into  classes,  with  the  class  boundaries  at  the  valleys 
of  the  histogram  (or  at  a break  in  slope,  leading  to  a shoulder).  Each 
peak  or  shoulder  would  then  be,  in  effect,  a modal  value  for  that 
spectral  band  for  a terrain  type.  Table  1 illustrates  the  product  of 
such  a process.  The  classes  were  chosen  subjectively;  others  may  choose 
somewhat  different  class  boundaries  in  some  cases. 

lU.  In  consideration  of  the  discussion  above,  a naturalistic  sub- 
division of  the  CCT  values  in  one  wavelength  band  is  defined  as  a "CCT 
value  class."  Thus,  a CCT  value  class  is  assumed  to  be  characterized 
by  a normal  distribution  of  CCT  values  around  a specified  modal  value, 
with  the  class  limits  some  fixed  standard  deviation  away  on  either  side 
of  the  modal  value. 

15.  Since  each  pixel  in  a Landsat  image  is  represented  by  four 
CCT  values  (i.e.  by  a CCT  value  set),  it  is  at  this  point  easy  to  leap 
to  the  conclusion  that  each  terrain  type  is  represented  by  a unique 
combination  of  CCT  value  classes.  However,  the  fact  that  not  all  wave- 
length bands  exhibit  the  same  number  of  CCT  value  classes  (e.g.,  band  !4 
has  7)  and  band  6 has  9)  illustrates  a critical  point;  namely,  not  all 
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of  the  CCT  values  in  a set  need  he  unique.  Only  the  set  as  a whole  is 
unique.  For  example,  one  terrain  type  might  be  characterized  by  values 
in  classes  1,  1,  1,  1 (bands  4,  5,  6,  and  T,  respectively),  and  another 
terrain  type  characterized  by  classes  1,  1,  1,  2.  Thus,  at  least  poten- 
tially, there  are  as  many  unique  spectral  "signatures"  as  permutations 
of  CCT  value  classes.  In  the  example  (Table  l),  there  may  be  4032  (7 
X 8 X 9 X 8 = 4032)  unique  combinations  of  classes. 

16.  For  future  discussion,  let  an  array  of  four  CCT  value  classes 
be  called  a "CCT  class  set."  Hence,  a CCT  class  set  is  assumed  to  be 
the  "spectral  signature"  of  a unique  combination  of  conditions  on  the 
ground. 

Relating  spectral 
signatures  to  terrain  types 

17.  The  central  questions  at  this  point  are:  Wiich  of  the  possi- 
ble combinations  are  "real"  (i.e.  originate  from  an  identifiable  terrain 
condition  or  type  on  the  ground),  and  how  may  the  relation  between  the 
spectral  signature  (i.e.  CCT  class  set)  and  the  terrain  type  be  estab- 
lished? The  following  two  basic  modes  by  which  these  questions  can 

be  answered  are  dependent  on  the  state  of  prior  knowledge : 

a.  An  "unguided"  mode  with  no  prior  knowledge  of  the  terrain 
type  spectral  signature  relation.  In  other  words,  the 
analyst  has  no  knowledge  of  the  actual  distributions  of 
terrain  types  on  the  groimd  or,  indeed,  of  what  types 
actually  exist  in  the  region.  Two  major  variants  of  the 
unguided  mode  are  known;  these  are  discussed  below. 

b.  A "guided"  mode  with  analyst  knowing  what  terrain  types 
exist  and  the  geographic  locations  of  some  presumably 
typical  or  characteristic  examples  of  each. 

18.  Ungulded  mode,  first  variant.  In  concept,  it  would  be  pos- 
sible to  construct  a map  showing  the  geographic  distribution  of  all 
pixels  characterized  by  any  specified  CCT  class  set.  For  example 
(Table  l),  all  class  1 CCT  value  ranges  (22-26,  band  4;  15-I8,  band  5; 
17-23,  band  6;  and  5-11,  band  7)  would  be  selected,  and  all  pixels  ex- 
hibiting any  set  of  values  encompassed  by  these  ranges  would  be  located 
in  geographic  space.  These  locations  would  then  (somehow! ) be  examined 
in  the  field  to  determine  what  terrain  condition  or  type,  if  any,  was 
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the  source.  This  process  would  be  repeated  for  each  of  the  possible 
permutations.  A situation  such  as  that  defined  by  the  Fif^ure  2 
histograms  would  mean  U032  such  exercises.  In  most  cases,  the  sheer 
magnitude  of  the  task  makes  this  approach  impractical. 

19-  There  is,  however,  a way  of  dramatically  reducing  the  size  of 
the  task.  It  depends  on  the  intuition  that  not  all  possible  permuta- 
tions of  CCT  value  classes  represent  actual  conditions  on  the  ground. 
Thus,  an  obvious  approach  is  to  determine  the  number  of  pixels  in  each 
permutation  of  CCT  value  classes.  Since  most  (if  not  all)  landscapes 
incorporate  only  a very  limited  number  of  distinct  terrain  types,  it 
would  be  expected  that  the  result  would  be  a relatively  small  number  of 
CCT  class  sets  having  a very  large  number  of  pixels.  It  would  also  be 
expected  that  a relatively  large  number  of  CCT  class  sets  would  show  a 
very  small  number  of  pixels.  The  latter  condition  would  be  expected  be- 
cause, purely  by  chance,  spurious  pixels  (see  paragraph  12c)  will  com- 
bine CCT  values  in  such  a way  as  to  meet  the  criteria  uf  one  or  another 
of  the  possible  CCT  class  sets. 

20.  At  this  point,  an  arbitrary  minimum  number  of  pixels  could  be 
selected  as  a cutoff  between  presumed  "real"  signatures  (those  deriving 
from  actual  terrain  types)  and  "false"  signatures  (those  deriving  from 
spurious  pixels).  When  this  procedure  is  actually  employed,  the  number 
of  CCT  class  sets  exhibiting  significantly  large  pixel  numbers  falls  to 
relatively  small  numbers,  normally  between  20  and  70.  This  is  still  a 
comparatively  large  number  to  field  check,  but  the  number  is  much  more 
tractable  than  the  original  several  thousand. 

21.  It  is  instructive  to  examine  the  reasons  for  the  persistence 
of  such  large  numbers  of  CCT  class  sets.  The  basic  reason  is  that  the 
sensor  ' sees"  and  records  only  radiation,  and  any  cliange  in  radiation 
Intensity  is  duly  recorded.  Therefore,  any  change  in  condition  that 
results  in  a change  in  reflectance  is  perceived  as  a new  spectral  signa- 
ture. On  the  other  hand,  human  perception  of  terrain  type  tends  to  be 
based  on  intrinsic  qualities  of  the  terrain.  A hardwood  forest  covering 
all  sides  of  a hill  is  perceived  as  "forest,"  uniform  over  the  entire 
hill,  however,  the  satellite  sees  the  slope  away  from  the  sun  quite 


differently  from  the  slope  toward  the  sun  and  thus  records  at  least  two 
(and  in  practice,  several)  signatures  from  the  same  hill.  The  conse- 
quence is  that  the  sensor  records  several  terrain  types  (i.e.  different 
CCT  class  sets)  for  the  feature  that  human  perception  and  experience 
insist  is  a single  homogeneous  terrain  type.  Since  it  is  the  human  needs 
that  must  be  satisfied,  the  problem  is  to  determine  which  of  the  sensor- 
perceived  signatures  are  derived  from  each  of  the  terrain  types  required 
by  the  human  analyst. 

22.  For  example,  suppose  the  analyst  wanted  to  map  the  distribu- 
tion of  only  three  terrain  types:  open  water,  forest,  and  nonforest. 

In  effect,  this  means  that  the  k032  potential  CCT  value  sets  (see  para- 
graph 15)  or  the  residual  20  to  JO  CCT  value  sets  (see  paragraph  20 ) 
must  be  allocated  among  only  those  three  categories.  In  the  absence  of 
any  information  on  either  the  location  of  actual  areas  exhibiting  the 
three  terrain  types,  or  information  on  the  probable  signatures  of  the 
three  types,  nothing  further  can  be  done.  Further  interpretation  re- 
quires some  process  by  which  the  relations  among  terrain  types  and  CCT 
value  sets  can  be  established. 

23.  As  described  above,  one  of  the  major  sources  of  multiple 
signatures  will  be  variations  in  reflectance  geometry.  This  fact  can  be 
exploited  to  assist  in  the  assembly  of  signatures  into  groups  or  arrays 
that  derive  from  a single  terrain  type  (as  perceived  by  the  human 
analyst).  When  a terrain  type  such  as  a forest  covers  an  entire  hill, 
the  CCT  value  sets  deriving  from  that  terrain  type  will  include  examples 
from  all  slope  angles  and  aspects  inherent  in  the  relations  among  topo- 
graphic configuration,  sun  zenith  angle  and  azimuth,  and  sensor  zenith 
angle  and  azimuth.  These  signatures  thus  form  a graded  sequence,  with 
the  constraint  that  the  same  reflectance  geometry  will  always  produce 
the  same  CCT  value  set.  A different  terrain  type  will  also  produce  a 
graded  sequence,  but  the  CCT  value  set  for  a given  reflectance  geometry 
will  be  different  from  that  deriving  from  the  first  terrain  type. 

2k,  In  theory,  the  procedure  would  involve  careful  placing  of  the 
geographic  locations  of  the  various  CCT  value  sets  on  a topographic  map. 
Any  specific  topographic  slope  degree  and  aspect  would  then  be  selected. 
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and  all  CCT  value  sets  recorded  from  each  occurrence  of  the  specified 
slope  would  be  noted.  If  there  are  only  two  distinct  terrain  types,  then 
there  will  be  two  CCT  value  sets  in  the  pixels  deriving  from  that  speci- 
fied slope.  If  there  are  three  terrain  types,  then  there  may  be  three 

CCT  value  sets,  but  only  if  all  three  terrain  types  occur  in  all  topo- 

graphic expressions.  In  any  event,  systematic  examination  will  result 
in  the  assembly  of  the  20  to  70  CCT  value  sets  into  a very  significantly 
smaller  number.  Each  aggregation  of  CCT  values  identified  in  this  manner 
can  be  regarded  as  a terrain  type,  and  at  this  point  a map  of  the  dis- 
tribution of  those  terrain  types  (which  are  actually  only  aggregations 
of  related  CCT  value  sets)  can  be  drawn.  A minimal  amount  of  field 
checking  will  now  establish  the  relation  between  ground  condition  (i.e. 
terrain  type)  and  aggregation  of  CCT  values. 

25.  Unguided  mode,  second  variant.  Let  it  be  assumed  that  the 

analyst  has  no  immediate  way  to  establish  relations  among  CCT  value  sets 

and  terrain  types  as  they  exist  on  the  ground.  Nevertheless,  there  are 

circumstances  in  which  the  analyst  would  like  to  map  the  distributions 
of  the  naturalistic  CCT  class  sets.  It  should  be  recalled  that  the  CCT 
class  sets  are  based  on  all  four  spectral  values,  hence  no  one  image 
(i.e.  no  one  Landsat  spectral  band)  can  be  used  to  properly  identify  a 
class  set.  Even  color  composite  images  are  not  particularly  helpful, 
since  the  human  eye  is  not  a reliable  spectrophotometer. 

26.  The  "unguided  mode,  second  variant"  procedure  begins  with  the 
selection  of  a number  of  sample  sites,  directly  from  a Landsat  image. 

The  sites  selected  would  normally  be  in  areas  of  homogeneity.  That  is, 
they  would  be  in  patches  in  which  the  gray  level  or  tone  of  the  image 
(i.e.  the  CCT  value)  was  reasonably  uniform  in  all  four  Landsat  bands 
for  at  least  several  pixels  in  all  directions.  This  minimizes  (but  does 
not  entirely  eliminate)  the  possibility  of  inadvertently  selecting  a 
spurious  pixel. 

27.  The  selected  pixels  would  then  be  separated  from  the  image 
array,  and  all  four  CCT  values  for  each  pixel  assembled  into  a CCT  value 
set.  Since  the  various  CCT  value  sets  are  all  assumed  to  be  members  of 
naturalistic  CCT  class  sets,  the  remaining  problem  is  to  determine  the 
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compositions  of  the  relevant  CCT  class  sets.  This  can  be  achieved  by 
any  one  of  several  forms  of  "cluster  analysis."  A rationale  for  such  a 
procedure  will  be  discussed  in  paragraphs  33-^3. 

28.  The  second  variant  is  constrained  by  a critical  problem:  sam- 
ple sites  must  represent  all  of  the  naturalistic  classes  of  CCT  values. 

If  they  do  not,  then  many  pixels  that  exhibit  CCT  value  sets  will  not 
belong  to  any  of  the  recognized  terrain  types  and  so  cannot  be  classi- 
fied. Thus,  it  seems  better  to  err  in  the  direction  of  abundance  rather 
than  to  select  too  few  and  have  residual  unclassified  areas. 

29.  Guided  mode.  The  guided  mode  (i.e.  in  which  the  analyst  has 
at  least  some  prior  knowledge  of  the  terrain  types  in  the  region  as 
well  as  knowledge  of  the  locations  of  at  least  some  characteristic 
examples  of  each)  is  a much  more  direct  method  of  exploiting  the  natural- 
istic CCT  value  classes.  This  mode  depends  on  the  fact  that  the  analyst 
can  establish  a direct  relation  among  terrain  types  and  CCT  value  sets 

or  CCT  class  sets.  In  principle,  a pixel  falling  inside  the  boundaries 
of  a known  example  of  a terrain  type  will  exhibit  a CCT  value  set  that 
will  be  characteristic  of  one  of  the  naturalistic  CCT  class  sets,  as 
defined  by  histograms  of  the  type  shown  in  Figure  2.  For  example,  a CCT 
value  set  of  25,  195  39,  30  (bands  k,  5,  6,  and  7,  respectively)  would 
fall  into  CCT  value  classes  (Table  l)  1,  2,  6,  and  6.  If  the  class 
boundaries  actually  represent  realistic  natural  divisions  among  terrain 
types,  the  CCT  value  set  is  a member  of  a group  of  CCT  value  sets,  each 
value  of  which  may  vary  within  the  limits  of  the  stated  classes.  That 
is,  the  value  of  band  h may  vary  between  22  and  26,  the  value  of  band  5 
between  I8  and  21,  etc. 

30.  When  this  has  been  done,  several  CCT  value  sets  will  be  dis- 
covered that  are  actually  members  of  the  same  population  (or  CCT  class 
set),  as  defined  by  the  naturalistic  classes  determined  from  the  histo- 
grams. Furthermore,  since  the  sample  sites  were  selected  because  they 
represented  the  various  terrain  types  of  interest  to  the  analyst,  the 
naturalistic  CCT  class  sets  associated  with  each  terrain  type  can  be 
readily  identified.  For  example,  variations  in  reflectance  geometry 
might  be  found  to  result  in  five  or  six  naturalistic  CCT  class  sets 
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being  representative  of  the  same  desired  terrain  type.  These  CCT  class 
sets  can  now  be  used  as  numerical  criteria  by  which  all  other  pixels  in 
the  scene  can  be  quickly  evaluated.  In  theory,  every  pixel  in  the 
scene,  the  CCT  values  of  which  fall  within  the  stated  limits,  will  be 
evidence  of  the  presence  of  the  terrain  type  on  the  ground.  With  the 
identification  of  all  pixels  meeting  the  specified  criteria,  a map  show- 
ing their  distribution  can  be  drawn;  and  that  map  will,  in  effect,  be  a 
representation  of  the  areal  distribution  of  the  specified  terrain  type. 

31.  A variant  of  this  method  assumes  the  existence  of  naturalistic 
CCT  class  sets  but  does  not  actually  use  the  histograms.  Instead,  it 
assuiries  that  the  CCT  value  exhibited  by  a sample  pixel  is  the  modal 
value  of  a set  of  values  that  is  normally  distributed  about  that  value. 

An  examination  of  the  sample  histograms  (Figure  2)  and  the  naturalistic 
classes  that  can  be  derived  from  them  reveals  that  the  shapes  of  the 
distributions  could,  in  most  instances,  be  generated  by  normal  distri- 
butions only  four  to  six  CCT  values  wide  (i.e.  two  or  three  values  on 
either  side  of  the  "peak").  Therefore,  it  is  not  entirely  unreasonable 
to  assume  that  the  value  of  the  sample  point  represents  the  modal  value 
of  a normal  distribution  with  a standard  deviation  of  1.5"  The  values 
away  from  the  mode  represent  system  errors,  as  described  in  paragraph  12. 

32.  In  effect,  each  CCT  value  set  is  assumed  to  define  the  modal 
values  of  a CCT  class  set.  However,  if  for  no  other  reasons  than  the 
system  errors,  CCT  value  sets  from  different  sample  areas  within  the 
same  terrain  type  may  exhibit  slightly  different  CCT  value  sets.  These 
closely  related  sets  shoixld  clearly  be  assembled  into  an  aggregation 
that  is,  collectively,  a CCT  class  set.  This  involves  another  aspect  of 
cJ  istex'  analysis  that  is  discussed  in  the  following  paragraphs. 

Cuu^ .. . ’ ..  ysis 

33.  To  be  of  maximum  utility  in  meeting  the  objective  stated  in 
paragraph  3,  the  cluster  analysis  procedure  must  be  simple,  rapid,  and 
under  operator  control  during  all  phases.  The  primary  reason  for  the 
last  requirement  is  that,  in  at  least  one  context  (see  paragraphs  29-32), 
the  analyst  will  be  trying  to  "force"  the  CCT  value  sets  into  a very 
limited  number  of  quite  arbitrary  groupings.  For  example,  a region 
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might  actually  display  20  or  more  "real"  GOT  class  sets,  but  the  analyst 
might  desire  that  his  final  product  consist  of  a map  with  only  three  or 
four  categories.  It  is  obvious,  in  such  a case,  that  each  mapped  cate- 
gory will  contain  several  of  the  naturalistic  CCT  class  sets. 

3U.  In  practice,  a CCT  value  set  obtained  for  a pixel  that  was 
selected  as  representing  a particular  terrain  type  cannot  be  assumed  to 
consist  of  CCT  values  representing  the  modal  values  of  the  CCT  class  set 
to  which  it  belongs,  since  instrument  errors  and  other  causes  can  alter 
the  theoretical  values  to  some  degree.  Further,  the  degree  of  variation 
cannot  be  known  & priori , and  thus  the  assumption  of  a fixed  standard 
deviation  (see  paragraph  31)  cannot  be  relied  upon.  Indeed,  there  is  no 
reason  to  assume  that  all  terrain  types  will  exJiibit  the  same  amount  of 
variability  in  their  reflectance  characteristics,  and  thus  it  may  be 
confidently  expected  that  some  naturalistic  CCT  value  classes  will  be 

broader  than  others.  An  examination  of  the  subjectively  selected  class  ' 

ranges  identified  in  Table  1 supports  this  exi)ectation ; some  classes  are 

only  two  CCT  values  wide,  and  others  ar?  as  much  as  nine.  y 

35.  Objective  criteria  could,  in  principle,  be  obtained  by  se- 
lecting many  pixels  representing  the  seune  terrain  type  ajid  considering 
them  to  be  members  of  the  s!ime  population,  with  values  arranged  in  a 
Gaussian  distribution.  The  simplest  method  of  achieving  this  would  aj - 
parently  be  to  select  a sample  pixel  representing  a ktu  wn  terrain  type 
and  then  to  enlarge  the  sample  to  include  mi  ai'ray  pixels  neigliboring 
the  sample  pixel.  If  only  the  immediate  neight'oi  s are  included,  llie 

combined  sample  would  consist  of  a 3-  by  i-pixel  array,  whicii  is  enougli  | 

to  provide  a sufficient  estimate  of  mode  and  standard  devia;  i- n 1 1’  tlie  I 

population  is  internally  homogeneous,  lliut  is,  tiie  assuni}  tion  is  valid 
only  if  the  va^-iability  results  from  sucii  sources  as-  inst  rument  cwors , 
small  and  randomized  variations  in  topiographic  slope,  etc.  If  ttie  sam- 
ple areas  are  chosen  with  care,  tlie  assumption  is  likely  to  be 
acceptable. 

36.  However,  it  is  not  necessary  to  proceed  entirely  on  faith. 

If  the  sample  3-  by  3-pixel  array  is  inadvertently  positioned  such  tliat 
some  of  the  pixels  are  in  one  terrain  tyjie  and  the  remainder  in  anottier. 
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the  CCT  values  of  at  least  one  of  the  spectral  bands  will  show  either 
a bimodal  distribution  or  a strongly  skewed  distribution.  Thus,  the 
obvious  step  would  be  to  subject  the  CCT  values  from  each  wavelength 
band  to  a normalcy  analysis.  If  the  distributions  of  values  in  all 
wavelength  bands  prove  to  be  essentially  Gaussian  (i.e.  normal),  then 
it  may  be  safely  assumed  that  all  nine  pixels  represent  a single  terrain 
type. 

37-  If  several  sample  3-  by  3-pixel  arrays  were  selected  in  the 
same  terrain  type,  ideal  theory  would  dictate  that  the  CCT  class  sets, 
as  defined  by  the  standard  deviations  of  the  3-  by  3-pixel  samples, 
should  overlap  perfectly.  However,  this  cannot  be  expected  in  practice, 
because  of  small  local  variations  in  ground  conditions , random  instru- 
ment error,  and  the  like.  In  view  of  the  practical  realities,  simi- 
larity was  defined  to  exist  between  two  CCT  class  sets  if  all  CCT  value 
classes  in  the  two  CCT  class  sets  overlap.  For  example,  suppose  that 
three  3-  by  3-pixel  arrays  have  been  selected  in  a region  and  that  the 
Landsat  band-U  mean  and  standard  deviations  and  CCT  class  set  limits 
are  as  follows: 


Site  No. 

Mean 

Standard  Deviation 

CCT  Class 
Set  Limits 

1 

25.0 

0.8 

2lt.2-25.8 

2 

23. h 

0.9 

21.5-26.3 

3 

32.0 

l.T 

30.3-33.7 

On  the  basis  of  the  criteria  for  similarity  previously  defined,  sites 
1 and  2 are  similar,  because  their  CCT  class  sets  overlap;  but  site  3 is 
not  similar  and,  therefore,  is  presiunably  a member  of  a different  popu- 
lation. That  is,  it  is  presumably  representative  of  a terrain  type 
other  than  that  designated  by  sites  1 and  2. 

38.  It  should  be  noted  that  the  establishment  of  similarity 
between  two  sites  may  result  in  a broadening  of  the  effective  CCT  class 
set.  For  exiunple,  if  only  site  1 in  the  tabulation  in  the  previous 
paragraph,  is  considered,  the  CCT  class  set  limits  for  tiie  terrain  type 
represented  by  site  1 are  21+.2-25.8.  However,  since  tiie  CCT  class  set 
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limits  of  site  2 overlap,  the  two  sites  are  defined  as  denoting  the  same 
terrain  type.  This  means  that  the  CCT  class  set  limits  representing  that 
terrain  type  have  been  broadened  to  include  both  sites  1 and  2 and  are, 
therefore,  2^*. 2-26. 3. 

39-  The  implications  of  this  are  important.  First,  consider  a 
natural  slope  smoothly  curved  such  that  a part  of  it  is  at  right  angles 
to  incident  solar  rays  and  thus  is  illuminated  at  maximiom  intensity, 
while  other  parts  are  at  progressively  greater  angles  until  the  slope  is 
parallel  to  the  incident  rays.  The  result  is  that  the  reflected  radia- 
tion will  exhibit  a smooth  continuum  of  values  from  maximiun  to  minimum 
in  all  wavelength  bands.  However,  the  radiance  values  in  each  band  will, 
for  all  practical  purposes,  retain  the  same  relative  ratios  with  respect 
to  each  other.  A pixel  grid  placed  over  the  slope  will  result  in  a 
series  of  CCT  value  sets  that  are  generically  related.  If  these  CCT 
value  sets  are  treated  in  the  manner  described  in  paragraphs  36-38,  the 
product  is  likely  to  be  a series  of  CCT  class  sets  that  sequentially 
overlap  and  fom  a related  series,  all  of  which  represent  the  same  ter- 
rain type.  'Thus,  to  some  degree,  the  vagaries  of  reflectance  geometry 
are  compensated  for  by  the  method  of  defining  similarity. 

The  same  situation  will  be  true  for  other  terrain  types,  of 
course,  but  it  must  be  recalled  that  it  is  the  relative  ratios  of 
radiance  values  that  determine  the  uniqueness  of  each  terrain  type, 
rnis  attribute  should  allow  a number  of  terrain  types  to  be  differ- 
entiated, more  or  less  independently  of  reflectance  geometry. 

^1.  The  second  important  implication  is  that  it  is  entirely  pos- 
sible to  have  three  CCT  class  sets,  in  which  set  1 is  related  to  set  2 
and  set  2 to  set  3,  but  in  which  an  independent  comparison  of  sets  1 and 
3 would  not  indicate  a similarity.  For  example,  suppose  the  analyst 
selected  four  sample  sites  in  a terrain  type  and  obtained  CCT  class  sets 
for  each  (see  paragraphs  36-38).  Figiu’e  3 illustrates  a possible  com- 
bination. If  a "similarity"  matrix  is  assembled  (Figure  3),  it  will  be 
noted  that  CCT  class  set  1 is  similar  to  sets  2,  3,  and  h,  but  that  set 
2 is  not  similar  to  set  3.  Tlie  reason  is  that  the  CCT  class  sets  in 
band  6 do  not  overlap  in  sets  2 and  3.  Despite  the  fact  that  sets  2 and 
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CCT  Class  Sets 


Set 


No . 
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Band  5 
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1 
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30.0-32.8 

2 

2U. 8-26.7 

18.6-20.1 

39.O-i4O.l4 

30.0-33.0 

3 

23.6-25.0 

17.8-19.2 

36.8-38.8 

29. 2-31.1* 

k 

2l*.  5-26.7 

18.1-19.8 

37.2-39.6 

30.1-32.7 
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123 

Matrix 

1. 
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" CCT  Value 
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X 

X 

X 

X 

Set 
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2 

X 

X 

X 

Band  1* 

21*.  2-26. 7 

3 

X 

X 

X 

Band  5 

Band  6 

18.1-20.3 

36.8-1*0.1* 

1* 

X 

X 

X 

X 

Band  7 

29.2-33.0 

Note:  X indicates 

similarity 


Figure  3.  Stretched  CCT  value  ranges  for 
a similarity  matrix 

3 are  not  similar  to  each  other  on  a one-to-one  basis,  they  are  defined 
in  this  study  as  similar  by  virtue  of  their  mutual  relationship  to  the 
entire  similarity  matrix.  To  express  analytically  the  type  of  similar- 
ity shared  by  all  four  sample  sites,  the  CCT  value  rtmges  of  each  site 
must  be  extended  or  stretched  in  each  band  to  the  same  limits.  The  re- 
sulting stretched  ranges  (Figure  3)  then  represent  the  CCT  class  set  of 
the  similarity  matrix  as  a whole. 

h2.  Finally,  it  will  be  recalled  that  the  limits  of  the  CCT  class 
sets  were  defined  arbitrarily  as  being  one  standard  deviation  away  from 
the  mean  on  either  side.  Tliere  is,  of  course,  no  reason  in  physics  why 
that  value  should  possess  any  intrinsic  merit.  If  the  class  ranges  (or 
boundaries)  are  set  farther  from  the  mean,  the  result  will  obviously  be 
that  more  CCT  class  sets  will  meet  the  similarity  criteria.  Thus,  ad- 
justments of  the  position  of  the  boimdaries  of  the  CCT  class  sets  will 
result  in  various  degrees  of  "tuning"  of  the  similarity  categories. 

Practice 


1*3.  'fhe  following  discussion  of  the  semiautomated  procedure  is 
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based  on  the  section  on  theory  (see  paragraphs  8-42)  and  on  ground  truth 
data  obtained  in  an  11-  by  l4-km  study  area  centered  approximately 
4o  km  northeast  of  Vicksburg,  Mississippi  (Figure  4).  The  study  area 
was  selected  to  include  a representative  section  of  the  loess  hills, 
forming  the  eastern  wall  of  the  Mississippi  floodplain,  and  a section  of 
the  floodplain,  including  an  oxbow  lake  and  a number  of  other  floodplain 
features.  Figure  5 is  a mosaic  made  of  conventional  panchromatic  air 
photos  of  this  study  area.  The  soil  type  categories  were  transferred  to 
the  photograph  from  the  U.  S.  Department  of  Agriculture  (USDA)  Soil  Con- 
servation Service  Soil  Survey  of  Yazoo  County^  (see  paragraph  50). 

44.  Plate  1 is  a generalized  flow  diagram  of  the  steps  involved 
in  the  procedure  described  below  (see  paragraphs  62-107  for  a detailed 
discussion) . 

45.  Since  the  study  area  undergoes  profoixnd  seasonal  changes  and 
a general  procedure  that  would  be  workable  in  any  season  was  sought,  the 
intent  was  to  use  one  data  set  (specifically  the  13  October  1975  data) 
to  develop  the  procedure  and  then  to  apply  that  procedure  rigorously  to 
the  study  area  as  represented  by  two  additional  Landsat  images , one  in 
summer  (ll  July  1974)  and  one  in  winter  (21  February  1975)*  The  July 

1974  imagery  data  were  obtained  with  Landsat  1,  and  the  February  and 
October  1975  images  with  Landsat  2.  All  three  scenes  were  carefully 
selected  to  be  free  of  both  haze  and  clouds,  thus  providing  optimiun  CCT 
values.  Table  2 gives  the  characteristics  of  the  imagery  data.  Note 
that  all  three  data  sets  were  of  approximately  the  same  overall  quality. 
The  only  obvious  differences  in  the  data  sets  are  in  terms  of  solar 
zenith  angle  and  azimuth.  Figure  6 is  an  image  formed  from  the  October 

1975  band-7  (O. 8-1.1  pm)  data,  indicating  the  location  of  the  study  area. 
Description  of  study  area 

46.  The  study  area  (Figure  5)  contained  several  different  terrain 
types,  including  agricultural  lands  (plowed  fields,  fields  with  crops  at 
different  stages  of  growth,  fallow  land,  etc.),  forested  land  (deciduous 
and  mixed  evergreen  trees),  grasslands  (summer  and  winter  pasture 
grasses,  fields  with  volunteer  grasses  and  other  weedy  plants),  water 
features  (rivers,  lakes,  canals,  wetlands),  and  some  low  topographic 
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Fig^e  6.  Image  provided  by  the  KROS  Data  Center  constructe  i 
irom  band-,  data  showing  location  ol  study  area; 

scene  22GI-155U0,  13  October  I975 
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areas  that  are  frequently  subjected  to  backwater  flooding.  It  covers 
small  sections  of  two  parts  of  the  Coastal  Plain  physiographic  province: 
(a)  the  Mississippi  Alluvial  Plain  and  (b)  the  East  Gulf  Coastal  Plain. ^ 
The  alluvial  soils  are  primarily  silts  (ML,  CL)  and  clays  (CH)  with  some 
sand;  whereas,  the  soils  in  the  East  Gulf  Coastal  Plain  are  primarily 
windblown  silts  (ML)  with  subordinate  amounts  of  very  fine  sand  and 
clay  (i.e.  loess  material). 

U7.  The  vegetation  in  the  area  consists  of  crops  of  soybeans, 
cotton,  and  bermuda  ajid  rye  pasture  grasses;  the  forested  areas  contain 
different  plant  densities  of  oak,  gum,  hickory,  poplar,  pine,  cedar, 
elm,  pecan,  and  cypress. 

i+8.  The  relative  effects  of  the  different  terrain  and  climatic 
factors  on  the  transfer  of  incident  radiant  energy  to  energy  that  has 
been  emitted  or  reflected  by  the  in  situ  terrain  materials  (soil,  rock, 
vegetation  material,  and  water)  and  man-made  features  (paved  areas, 
buildings,  etc.)  have  not  been  quantitatively  determined.  Nevertheless, 
there  are  some  basic  factors  that  almost  surely  act  together  in  some 
fashion  and  control  this  transfer  process.  These  factors  are  extremely 
critical  to  an  overall  procedure  that  will  allow  interpretation  of  ter- 
rain conditions  for  different  seasonal  and  climatic  conditions. 

U9.  Since  terrain  conditions  change  appreciably  with  time  and  are 
affected  significantly  by  local  climatic  conditions,  it  was  necessary  to 
quantitatively  describe  the  study  area  at  the  time  of  each  respective 
overflight.  Data  pertaining  to  soil,  topographic  and  flood  conditions, 
rainfall,  vegetation,  and  other  features  were  obtained  and  mapped  at  an 
appropriate  scale  for  comparison  with  the  conditions  being  interpreted 
from  the  Landsat  spectral  data.  The  mapped  data  were  portrayed  on  a 
photomosaic  at  a scale  of  1:1+000  (Figure  5)  that  was  made  up  from  1973 
air  photos  obtained  by  the  USDA. 

50.  Soil  conditions.  The  soils  within  the  study  area  were  mapped 
by  the  USDA  Soil  Conservation  Service^  and  consist  of  silty  loams  (ML), 
silty  clay  loams  (CL),  and  clays  (CH).  A soil  map  for  the  study  area  is 
superimposed  on  the  photomosaic  in  Figure  5*  The  soils  in  the  alluvial 
plain  (or  cropland)  are  composed  of  patches  of  ML,  CL,  and  CH;  whereas. 
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the  soils  in  the  upland  forested  areas  are  primarily  ML's. 

51.  Topographic  and  flood  conditions.  The  relief  within  the 
study  area,  as  portrayed  on  a U.  S.  Geological  Survey  (USGS)  1:62, 500- 
scale  map  with  5-ft  (1.52  m)  contours,  varies  from  about  85-115  ft 
(25.9-35.1  m)  above  mean  sea  level  (msl)  within  the  cropland  or  agri- 
culture region  to  about  115-130  ft  (35.1-39.6  m)  within  the  upland  hills 
that  occupy  the  remaining  portion  of  the  study  area. 

52.  The  study  area  borders  on  the  Yazoo  River,  a large  tributary 
of  the  Mississippi  River,  and  a portion  of  the  agricultural  area  is  sub- 
jected to  flooding  during  high  river  stages.  The  Corps  of  Engineers 
stage  records  for  the  Yazoo  River  gage  at  Satartia,  Mississippi,  in- 
dicate that  flooding  does  occur  within  adjacent  low-lying  areas  when 
the  stage  of  the  river  reaches  approximately  95  ft  (29  m)  msl.  The 
river  stage  data  from  the  Satartia  gage  for  the  period  10  January  197^ 
through  30  October  1975  (Figure  7)  indicate  that  flooding  most  likely 
occurred  within  parts  of  the  study  area  prior  to  and  diiring  two  of  the 
three  overflight  times  (ll  July  197^  and  21  February  1975). 

53.  Since  the  Yazoo  River  at  Satartia  did  attain  stages  that 
would  permit  flooding  (Figure  7)  in  low-lying  areas  during  two  of  the 
three  overflight  times,  it  was  anticipated  that  a detailed  topographic 
map  was  needed  in  the  event  that  it  became  necessary  to  delineate  low- 
elevation  areas  within  the  study  area.  The  best  available  map  was  the 
USGS  1:62, 500-scale  topographic  map.  The  topographic  map  was  constructed 
from  1963  aerial  photography  and  does  not  represent  the  small  changes 
and  alterations  that  have  been  made  to  the  area  to  make  the  low-lying 
areas  more  useful  for  agricultural  purposes.  Ditches  and  canals  have 
been  constructed  within  the  area  (Figure  8)  to  provide  for  drainage  dur- 
ing wet  periods  and  irrigation  during  dry  periods.  Figure  9 shows  the 
section  of  the  map  covering  the  study  area. 

Rainfall.  The  daily  rainfall  records  from  Germania,  Missis- 
sippi (site  in  the  middle  of  study  area.  Figure  9) > show  that  the  area 
received  some  rainfall  within  5-6  days  of  the  July  197 and  February 
1975  satellite  overflights,  and  no  rainfall  within  13  days  prior  to  the 
October  1975  cverflight  (Figure  10).  This  record  is  believed  to  be 
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Figure  8.  Canals  and  ditches  in  the  study  area,  8 May  1976 
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Figure  7.  Yazoo  River  stages  at  Satartia,  Mississippi 
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representative  of  the  rainfall  received  for  the  entire  study  area,  espe- 
cially when  compared  with  the  records  obtained  at  other  nearby  stations 
(Onward,  Rolling  Fork,  Vicksburg,  and  Yazoo  City).  Figure  11  presents 
two  years  (197^  and  1975)  of  monthly  rainfall  data  for  the  Germania 
station. 

55.  Classification  of  terrain  types.  The  in  situ  terrain  con- 
ditions for  the  three  overflight  periods  were  mapped  by  using  the  1973 
USDA  1:10,000  aerial  photography  and  ground  surveys.  Since  two  of  the 
overflights  (July  197^  and  February  1975)  occurred  before  the  time 
period  in  which  this  study  was  conducted  (September  1975-June  1976),  it 
was  impossible  to  obtain  some  of  the  ground  truth  data  that  would  have 
permitted  a more  complete  verification  of  the  Landsat  interpretation  and 
mapping  methodology  for  those  two  Landsat  data  sets;  however,  it  was 
possible  to  obtain  some  data  on  the  area  for  these  periods  pertaining  to 
crop  types  and  locations  by  contacting  various  landowners  within  the 
area. 

56.  The  terrain  conditions  within  the  study  area  were  mapped 
according  to  the  following  category  types: 
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Land-Use/Terrain  Types 


Agriculture:  foliated  cotton 
Agriculture:  defoliated  cotton 
Agriculture:  soybeans 

Agriculture:  oats 

Agriculture:  winter  wheat  (rye  grass) 

Agriculture:  volunteer  vegetation 
(no  planted  crops) 

Agriculture:  plowed  land  (no 
vegetation) 


F^  Floodplain  forest : 90  percent 

oak,  elm,  gum,  and  pecan;  10  percent 
cypress 
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hickory,  and  poplar;  15  percent 
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Q.  RAINFALL  IN  1974 


b.  RAINFALL  IN  1975 


Figure  11.  Monthly  rainfall  data  for  the  Germania  station  from 
January  197^  through  December  1975 
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Map  Symbol 


F 


3 


P 

P 

W 


1 

2 


Land-Use /Terrain  lYpes 

Upland  forest;  oak,  gum,  hickory, 
and  poplar 

Pastureland : bermuda  grass 
Pastureland:  bermuda  and  rye  grasses 
Water:  lakes,  rivers,  and  canals 


57-  The  maps  for  the  three  overflight  times  (Figures  12,  13,  and 
ill)  show  that  the  land-use  and  vegetation  types  changed  significantly 
from  the  first  overflight  (July  19T1+)  to  the  last  (October  1975)*  Not 
only  did  the  type  of  agricultural  crops  change,  but  changes  occurred  in 
some  of  the  floodplain  forested  areas.  Some  of  the  forested  areas  were 
cleared,  and  the  land  was  converted  to  agricultural  crops  and  pasture 
grasses  (Figure  15).  In  addition,  local  I’esidents  indicated  that  some 
of  the  study  area  had  been  flooded  at  the  time  of  the  July  1971+  and 
February  1975  overpasses. 

58.  Detailed  fieldwork  was  performed  at  the  time  of  the  October 
1975  overflight,  and  the  data  collected  were  used  to  prepare  the  terrain 
map  shown  in  Figure  lU.  Ground  photos  were  taken  of  most  of  the  sites 
at  the  same  time;  Figures  16-20  present  in  detail  some  of  the  signifi- 
cant terrain  types  that  occurred  within  the  study  area.  Figure  21  indi- 
cates the  locations  of  these  ground  photos. 

59.  To  obtain  additional  information  pertaining  to  the  terrain 
conditions  within  the  study  area,  low-altitude  air  photos  were  taken  on 
8 May  1976  (Figures  22-30)  and  '^re  located  as  indicated  in  Figure  31. 
These  air  photos  show,  in  detail,  the  complexity  of  the  terrain  condi- 
tions and  some  of  the  changes  that  had  been  made  to  the  landscape. 
Automated  data  station 

60.  An  automated  meteorological  station  (Figure  32)  was  in- 
stalled near  Germania  (coordinates  7.8-3.0,  Figure  9)  to  obtain  detailed 
information  on  incoming  solar  radiation,  rainfall,  wind  speed  and  direc- 
tion, and  air  and  soil  temperatures.  The  original  intention  had  been  to 
collect  such  data  for  at  least  a two-week  period  prior  to  the  October 
1975  Landsat  overpass.  However,  the  equipment  could  not  be  made  avail- 
able at  that  time,  and  thus  the  actual  period  of  record  extended  from 
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b.  Cropland 


Figure  15.  Areas  shovm  with  forest  removed  (since  air 
photos  dated  July  1973)  are  currently  used  (8  May  1976) 
as  pastureland  and  cropland  on  the  southeast  corner  and 
eastern  edge,  respectively,  of  the  large  forested  area 
on  east  side  of  Dump  Lake  (Figure  lU) 
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Close-up  of  leafy  surfaces 


Fifiure  17.  Mature  foliated  cotton,  terrain  type  A^,  13  October  1975 


Defoliated  cotton,  terrain  tyjio  A , IS  October  1075 
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Figure  19-  Recently  plowed  fields,  terrain  type  A_,  13  October  1975 


b . Uninowed 

Fif^ure  20.  Pastureland  berniuda  f^raar.,  terrain  type  1'  , 

13  October  197'.^ 
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Figure  24.  Study  area  as  viewed  south 
southeast  of  forested  area  surrounded 
by  cropland,  8 May  1976 


Figure  23.  Study  area  as  viewed  southwest  of  the  forested 
area  southeast  of  Dump  Lake,  8 May  1976 


Figure  25.  Study  area  as  viewed  east-southeast  of  Yazoo 
River  and  forest  that  lines  the  banks,  8 May  19T6 


area  as  viewed  northeast  of  forest  area; 
ide  of  Dump  Lake,  8 May  197^ 


Figure  27.  Study  area  as  vieved  northwest  of  Dump  Lake 
and  bordering  forest,  8 May  1976 


Figure  28.  Study  area  as  viewed  northeast 
of  complex  of  planted  fields,  pastureland, 
and  forest,  8 May  1976 


j|B 

Figure  30.  Study  area  as  viewed  north-northwest  of  small  open  areas 
(cropland)  surrounded  by  30-m  tall  deciduous  forest,  8 May  1976  at 
1700  hr;  note  effects  of  shadows 


Fij3ure  31 


Locations  and  directions  of  view  of  air  photo 
8 May  197o  (Figures  22-30) 


Figure  32.  Automated  meteorological  field  station 
as  installed  in  study  area  (coordinates  7- 8-3.0, 

Figure  5)?  28  October-18  November  1975 

28  October  to  l8  November  1975*  While  the  data  obtained  are  not  useful 
for  determining  such  conditions  as  soil  moisture  at  the  time  of  overpass 
they  were  nevertheless  helpful  as  a general  guide  to  meteorological 
characteristics  of  the  study  area  during  the  autumn. 

Basic  procedure  for 
guided  classification 

61.  The  actual  procedure  used  to  obtain  an  interpretation  of 
Landsat  imagery  data  was  an  interactive  process,  in  which  the  human 
analyst  performed  a number  of  nonquantif iable  operations  from  which 
quantitative  instructions  for  an  assortment  of  computer  programs  were  ob 
tained.  The  programs  then  very  rapidly  performed  a number  of  relatively 
simple  bookkeeping  operations  on  the  very  large  number  of  CCT  values 
comprising  a Landsat  image  data  set.  This  cycle  of  human  interpretation 
and  machine  operation  was  repeated  several  times.  The  general  scheme  is 
illustrated  as  a flow  diagram  (Plates  1 and  2).  Individual  segments  are 
shown  as  "blocks"  and  numbered  sequentially. 

62.  Block  1.  The  first  step  was  to  obtain  from  the  EROS  Data 


Center  (see  paragraph  2)  the  CCT's  for  each  of  the  three  Landsat  over- 
pass dates  selected  for  study  (Table  2).  The  data  for  each  overpass  con- 
sisted of  a pair  of  9-track  digital  magnetic  tapes  containing  multi- 
plexed CCT  values.  Two  tapes  together  comprised  a Landsat  scene  that 
covered  an  area  approximately  185  km  square. 

63.  Block  2.  The  study  area  selected  was  smaller  than  a complete 
scene;  therefore,  some  preprocessing  of  the  CCT  data  was  required.  The 
problem  was  to  extract  the  data  for  the  study  area  from  the  data  array 
covering  the  entire  scene.  In  principle,  it  should  have  been  possible 
to  simply  locate  the  study  area  carefully  on  a MSA-provided  1:1,000,000 
image  (one  of  the  user  product  options  available  from  the  EROS  Data 
Center)  and  then  the  corner  coordinates  by  measuring  in  from  the  edges 
of  the  image.  Unfortunately,  this  procedure  resulted  in  substantial 
error.  One  of  the  reasons  for  the  error  was  that  the  EROS-provided 
image  was  cropped  so  that  it  covered  somewhat  less  area  than  the  CCT 
data  array.  Since  the  amount  of  cropping  is  variable,  the  margins  of 
the  EROS-provided  images  cannot  be  used  to  precisely  fix  the  location  of 
a point  in  the  pixel  array. 

6h.  The  procedure  followed  was  to  write  an  image  of  the  scene, 
using  band-T  data  (because  of  the  clarity  with  which  water  surfaces  are 
evident),  and  simultaneously  making  approximate  corrections  for  earth 
rotation  effects  and  pixel  shape,  as  discussed  in  paragraphs  103-105  and 
in  Reference  3- 

65.  Blocks  3 and  U.  The  image  produced  after  making  these  ap- 
proximate corrections  is  sufficiently  close  to  the  geometries  of  con- 
ventional map  projections  that  it  was  easy  to  locate,  with  only  a small 
margin  of  error,  the  position  of  the  study  area.  Since  the  pixel  size 
used  to  produce  the  image  was  known  (50  pm  square),  and  the  edges  of  the 
image  represented  the  true  edges  of  the  CCT  data  array,  it  was  easy  to 
measure  from  the  image  margins  and  thus  approximately  fix  the  locations 
of  the  corner  coordinates  of  the  study  area  in  terms  of  pixel  and  scan 
line  numbers. 

66.  Block  5-  The  scan  line  arid  pixel  numbers  of  tiie  corner  co- 
ordinates were  used  as  the  basis  for  extracting  the  block  of  pixels 
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representing  the  study  area  from  the  complete  CCT  data  array.  The 
extracted  data  array  was  then  copied  onto  a new  tape,  which  was  called 
SA/CCT  (study  area  CCT). 


67.  Block  6.  The  next  step  in  the  procedure  was  to  write  an 
image,  at  an  enlarged  scale,  of  the  study  area,  using  selected  wave- 
length bands.  For  the  current  experiment,  band-T  data  were  chosen, 
chiefly  because  water  features  show  so  distinctly.  The  importance  of 
this  will  become  apparent  later  in  the  discussion.  In  the  current 
exercise,  conventional  air  photos  at  a scale  of  approximately  1:71,700 
were  used  for  certain  aspects  of  ground  truth,  and  so  the  images  were 
written  at  approximately  that  scale. 

68.  Block  7-  Since  a CCT  value  was  to  be  correlated  with  a point 

on  the  ground  (see  paragraphs  29  and  30),  the  location  of  the  point  (i.e.  | 

the  pixel  representing  that  point)  had  to  be  precisely  known.  An  error  j 

in  location  of  several  pixel  dimensions  (a  Landsat  pixel  is  a rectangle  ■ 

approximately  79.1  m long  in  a direction  parallel  to  the  orbital  path  of 

the  satellite  and  57-2  m wide  in  a direction  at  right  angles  to  the  I 

orbital  path)  could  have  resulted  in  an  attempt  to  correlate  a CCT  value  , 

set  with  the  wrong  terrain  type.  The  implication  of  this  was  that  the  i 

SA/CCT  had  to  be  corrected  to  the  geometry  of  an  air  photo  that  was  used  i 

as  the  basis  for  selecting  sample  items  (see  paragraph  29).  The  success  I 

of  the  procedure  for  achieving  geometric  accordance  depended  on  the 

existence  of  common  "transfer  points."  . 

69.  A critical  problem  was  the  selection  of  the  transfer  points.  ; 

They  had  to  be  such  that  they  could  be  located  in  the  Landsat  image  or  J 

data  set  regardless  of  the  season  at  which  the  data  were  obtained.  In  > 

the  current  instance,  data  for  three  widely  different  seasons  (summer,  i 

\ 

fall,  and  winter)  had  to  be  used.  It  must  be  recalled  that  the  transfer  j 

I 

points  must  also  be  clearly  recognizable  on  the  air  photo  as  well  as  in  : 

the  Landsat  data.  The  following  features  usually  provide  good  transfer  i 

i 

transfer  points:  J 

a.  A small  (approximately  2.5-ha)  roughly  circular  lake  | 

surrounded  by  dry  land  or  contiguous  vegetation.  j 

b.  Intersection  of  two  large  (i.e.  50-m  wide)  streams. 
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c_.  The  Junction  of  a large  stream  and  a lake. 

Intersection  of  two  wide  roads  or  airfield  runways. 

70.  The  large-scale  image  of  the  study  area  (Block  6)  and  an  air 
photo  of  approximately  the  sfime  scale  were  studied  carefully.  Then 
three  transfer  points,  A,  B,  and  C,  were  selected  (P'igure  33).  They  are 
not  ideally  located;  it  would  have  been  better  had  they  been  more  widely 
separated.  However,  as  it  happens,  the  southern  end  of  the  study  area 
offered  few  opportunities  for  suitable  transfer  points.  Low-altitude 
oblique  air  photos  (Figure  3^)  illustrate  the  kinds  of  features  at 
transfer  points  A and  C.  The  features  were  identifiable  in  the  Lar.dsat 
data  sets  for  July  197^  as  well  as  February  and  October  1975.  Fig- 
ure 35  presents  a 10-  by  10-pixel  matrix  of  the  band-7  OCT  values  around 
transfer  point  A.  Note  that  in  each  data  set  the  junction  of  the  stream 
with  the  lake  bridge  shows  as  a higher  value  than  the  "lake"  pixels  and 
that  the  small  wide  spot  in  the  stream,  just  before  it  enters  the  lake 
and  south  of  the  junction,  shows  as  a value  significantly  lower  than  the 
surrounding  dry  land.  The  junction  is,  of  coirrse,  readily  identifiable 
on  the  air  photo. 

71.  Block  8.  The  problem  now  was  to  stretch  (or  "rubber  sheet") 
the  GOT  pixel  array  in  such  a way  that  pixels  representing  transfer 
points  are  located  at  the  same  coordinates  as  their  correlative  points 
on  the  air  photo.  The  procedure  was  to  assign  coordinate  positions  to 
the  air-photo  locations  of  the  transfer  points  by  first  designating  one 
of  the  transfer  points  as  the  origin  of  a convenient  but  arbitrary 
Cartesian  coordinate  system  and  then  determining  the  X and  Y coordi- 
nates of  the  other  points  with  respect  to  this  origin.  Likewise,  the 
transfer  points  within  each  Landsat  data  set  were  located  by  determining 
the  number  of  pixels  (N  ) in  the  X direction  and  (N  ) in  tlie  Y 

A I 

direction  of  each  transfer  point.  The  map  coordinates  of  the  transfer 
points  were  then  related  to  the  pixel  coordinates  within  each  Landsat 
data  set  by  the  following  two  equations: 


X. 

1 

cos 

®f  - 

Y. 

1 

sin  0^  = N^. 

°Xi 

X. 

1 

sin 

«f " 

Y. 

1 

cos  0^  = 

^i 
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SAMPLE  SITES 


APPROXIMATE  SI7E  OF  S-HY  3-PIXEI 
ARRAYC.Q')  is  34  BY  2 4 MM. 


Location  of  transfon  points 
( + ) in  study  art'a 


0.  11  JULY  1974 


b.  21  FEBRUARY  1975 


c.  13  OCTOBER  1975  d.  AIRPHOTO 

LEGEND 

© TRANSFER  POINT  A 

P’if'ure  35.  B.'ind-T  CCT  value;;  around  transfer  point  A (all  three 
overpar.se;; ) with  air  jhioto  for  oompari;;on 


where 


X.  and  Y.  = the  X and  Y air-photo  coordinates  of  each  respec- 
tive transfer  point 

0^  = Landsat  flight  path  angle  with  respect  to  true  north 

= number  of  X-coordinate  pixels  of  the  transfer  point 
within  the  Landsat  data  set 

= width  of  the  Landsat  pixel 

= number  of  Y-coordinate  pixels  of  the  transfer  point 
within  the  Landsat  data  set 

= length  of  the  Landsat  pixel 

Equations  representing  each  respective  transfer  point  were  written  and 
then  simultaneously  solved  for  the  flight  path  angle  and  the  size  (D 

X 

and  D^)  of  the  pixel  within  each  Landsat  data  set.  The  values  of  0^ 
and  pixel  size  for  the  three  Landsat  data  sets  are  listed  in  Table  3. 

The  average  values  for  were  then  computed  and  used 

in  the  equations  above  to  bring  the  terrain  reference  system  and  CCT 
value  data  for  each  data  set  into  geometric  accordance. 

72.  The  number  of  transfer  points  needed  to  bring  Landsat  data 

sets  into  geometric  accordance  depends  on  the  size  of  the  area  being 

considered  and  the  relief  within  the  area.  It  is  estimated  that  in  most 

regions  in  the  United  States  between  12  and  20  points  would  probably  be 

2 

required  for  an  area  the  size  of  a NASA  scene  (“3^,000  km  ).  Reference  8 
describes  a procedure  for  bringing  Landsat  data  into  geometric  accordance 
using  up  to  20  transfer  points. 

73.  Block  9.  The  next  problem  was  to  select  the  sample  sites 
(see  paragraph  29)  to  be  used  in  the  "guided"  interpretation  mode.  How- 
ever, before  making  the  selection,  at  least  a few  examples  of  each 
terrain  type  had  to  be  actually  located.  In  the  present  exercise,  this 
was  accomplished  by  a combination  of  air-photo  interpretation,  using  the 
photomosaic  in  Figure  5,  and  ground  truth  data,  as  described  in  para- 
graphs U6-60.  The  initial  classification  was  comprised  of  only  three 
terrain  types;  open-water  surfaces,  wooded  land  (or  forest),  and  non- 
wocded  (or  open)  land.  On  the  basis  of  all  available  information,  30 
sample  sites  were  selected  (Figure  33)  that,  it  was  hoped,  would  include 
every  significant  variation  of  condition  in  the  study  area.  As  noted  in 
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Figiire  33,  a Landsat  pixel  is  an  area  about  1.1  mm  long  and  0.8  mm  wide, 
with  the  long  dimension  north-south.  Thus,  a 3-  by  3-pixel  array  (see 
paragraph  36)  is  a rectangle  3.^  mm  N-S  and  2.h  mm  E-W.  Each  sample 
site  was  selected  at  the  center  of  an  area  of  homogeneity  as  perceived 
on  the  air  photo  and  in  an  area  as  large  as  a rectangle  3.^  by  2.k  mm 
(Figure  33).  By  the  use  of  available  ground  truth  data  and  air-photo 
interpretation,  the  terrain  conditions  at  each  site  during  each  satellite 
overpass  were  established  (Table  h). 

T4.  Blocks  10  and  11.  A 3-  by  3-pixel  array  of  values  around  each 
sample  site  was  isolated  from  the  SA/CCT  and  subjected  to  a "normalcy" 
analysis  to  determine  actual  homogeneity.  Table  5 gives  the  mean  and 
standard  deviations  of  the  COT  values  for  each  wavelength  band  for  the 
October  1975  overflight.  As  noted,  the  standard  deviations  vary  consid- 
erably; for  example,  the  standard  deviations  are  quite  large  (greater 
than  2.0)  for  sites  17  and  29.  Also,  Figure  33  shows  that  there  are 
situations  such  that  a 3-  by  3-pixel  array  cannot  quite  be  fitted  into 
an  area  of  complete  homogeneity.  These  sites  would,  of  coxu'se,  be 
expected  to  exhibit  more  variance  than  others.  Ideally,  the  variance 
exhibited  by  a 3-  by  3-plxel  array  should  be  no  larger  than  the  variance 
produced  by  normal  instrument  error.  Since  the  nomal  instrument  error 
is  about  +2  percent,  the  ideal  standard  deviation  would  be  +1.26  CCT 
units  (maximum  CCT  value  of  63  multiplied  by  0.02,  the  normal  Landsat 
instrument  error).  The  assumption  is  that  any  3-  by  3-pixel  array 
exhibiting  a standard  deviation  of  +^1.26  would  be  likely  to  represent  an 
absolutely  homogeneous  terrain  type.  In  practice,  such  homogeneity  can- 
not be  expected  because  of  small-scale  randomized  variations  in  surface 
geometry,  soil  moisture,  etc.  Since  there  is  no  way  of  telling  ahead 
of  time  how  large  such  variations  are  likely  to  be,  the  analyst  should 
be  permitted  to  control  the  amount  of  variance  that  is  used  to  Judge  the 
homogeneity  of  a sample  site,  with  the  proviso  that  the  variance  should 
in  no  case  be  less  than  about  jil.3  (+1.26  rounded  to  the  nearest  0.1 ). 
This  controlled  parameter,  which  was  used  to  determine  the  acceptability 
of  the  3-  by  3-pixel  arrays,  will  henceforth  be  called  the  "screening 
window."  Because  some  degr-’e  of  site  variance  was  anticipated,  ii.5  CCf 
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units  were  initially  accepted  as  the  screening  window.  Those  sample 
sites  that  exhibited  standard  deviations  of  greater  than  +1.5  CCT  units 
in  one  or  more  of  the  four  spectral  bands  were  eliminated  from  further 
consideration.  The  assumption  behind  this  decision  was  that  a larger 
standard  deviation  would  be  prima  facie  evidence  that  two  or  more  terrain 
types  were  represented  in  the  sample.  This  decision  drastically  reduced 
the  number  of  samples  used  in  analysis;  for  example,  only  17  sites  re- 
mained out  of  the  original  30  for  the  October  1975  data. 

75.  The  specification  of  a particular  screening  window  as  the 
criterion  for  homogeneity  is,  of  course,  a specification  of  the  amount 
of  variance  that  will  be  permitted  in  a sample  site.  If  the  acceptable 
screening  window  is  made  larger,  more  sample  sites  will  be  found  to  be 
acceptable.  Table  6 illustrates  the  effect  of  increasing  the  screening 
window  on  the  number  of  accepted  sample  sites. 

76.  Examination  of  the  photomosaic  (Figure  3.3)  sheds  little  light 
on  the  reasons  for  the  variances  or  lack  thereof,  since  some  sites  that 
appear  on  the  photomosaic  to  be  very  homogeneous  actually  exhibit  quite 
large  variances.  It  seems  that  subjective  analysis  of  air  photos  will 
not  guarantee  the  selection  of  homogeneous  sample  sites.  Therefore,  the 
best  approach  would  be  to  select  a large  number,  expecting  that  many 
will  prove  unsuitable,  at  least  initially. 

77.  Block  12.  The  initial  terrain  classification  was  to  consist 
of  three  categories:  water,  forest,  and  nonforest  (see  paragraph  73). 
Thus,  each  of  the  sample  sites  that  survived  the  normalcy  analysis  (i.e. 
that  met  the  screening  window  criterion  (Blocks  12.1  and  12.2,  Plate  2)) 
must  be  assumed  to  be  a member  of  a "cluster"  of  sample  sites,  each  of 
which  represents,  in  some  degree,  one  of  the  three  terrain  types.  Thus, 
the  first  truly  subjective  decision  must  rest  on  the  apparent  represen- 
tativeness of  the  selected  sample  sites  (Block  12.3,  Plate  2).  If  only 
two  or  three  survive  the  normalcy  screening,  and  the  terrain  is  known 

to  be  complex  in  any  sense,  then  it  may  be  that  the  number  of  selected 
sample  sites  is  too  small.  In  this  event,  the  obvious  step  is  to  en- 
large the  screening  window  so  that  additional  sample  points  are  included 
( Table  6 ) . 
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78.  For  the  October  overpass  data,  it  was  decided  that  the  se- 
lected sites  adequately  represented  the  desired  terrain  types.  The  next 
step  was  to  test  the  selected  samples  for  similarity.  Ideally,  all  of 
the  samples  denoting  one  terrain  type  should  exhibit  a relatively  close 
degree  of  radiometric  similarity.  In  practice,  this  means  that  the  CCT 
values  of  all  of  the  samples  designating  a desired  terrain  type  should 
show  a sort  of  family  resemblance.  Furthermore,  the  CCT  values  repre- 
senting the  sEunples  in  a "terrain"  type  should  tend  to  "cluster"  about  a 
common  value,  and  the  cluster  for  each  terrain  type  should  be  discrete, 
i.e.  should  not  overlap. 

79*  The  first  step  in  the  .esting  process  was  to  select  a "vari- 
ance" window  (Block  12.4,  Plate  2).  Ideally,  the  smallest  variance  win- 
dow should  be  equal  to  or  larger  than  the  maximum  instrument  error,  and 
thus,  when  dealing  with  Landsat  data,  it  should  be  equal  to  or  greater 
than  ^1.26  CCT  units.  However,  the  selection  of  the  variance  window  is 
entirely  subjective  and  must  be  made  by  the  analyst.  Generally,  it 
appears  to  be  better  to  first  choose  an  initial  value  relatively  close 
to  the  limiting  minimum.  This  variance  window  (which  is  both  positive 
and  negative)  is  then  algebraically  added  to  the  mean  value  of  all  four 
spectral  bands  of  the  selected  sites  (Block  12.5,  Plate  2).  The  result- 
ing values  are  then  rounded  to  the  nearest  integer,  simply  because  the 
CCT  values  are  always  integers.  Table  7 presents  the  effects  of  this 
process,  using  three  different  variance  windows.  The  two  integer  values 
that  bracket  each  mean  value  represent  class  limits  for  each  spectral 
band.  In  effect,  the  concept  is  that  all  CCT  values  that  fall  between 
or  on  the  limiting  values  are  assumed  to  be  equivalent.  For  example, 

CCT  values  of  15,  l6,  or  17  in  spectral  band  4 (Table  7)  would  be  as- 
sumed to  be  members  of  a population  that  included  site  1,  given  a vari- 
ance window  ^1.3  CCT  units.  However,  note  that  increasing  the  size  of 
the  variance  window  tends  to  broaden  the  width  of  the  class.  In  physical 
terms,  this  means  that  a wider  range  of  CCT  values  will  meet  the  limiting 
criteria  and  are  thus  categorized  as  belonging  to  the  class. 

80.  The  next  procedural  step  (Block  12.6,  Plate  2)  was  to  con- 
struct a "cluster  matrix"  (Figure  36)  for  the  13  October  data,  using  as 


SITE  NUMBERS 
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Figure  36.  Cluster  matrix,  13  October  1975 
data,  with  variance  window  of  ;^1.5  OCT  units 

a first  choice  a screening  window  and  a variance  window  of  jd  . 5 CCT 
units.  The  matrix  is  simply  an  orthogonal  array  with  the  identification 
numbers  of  the  selected  sites  arranged  in  numerical  sequence  along  both 
abscissa  and  ordinate. 

8l.  To  "fill"  the  matrix,  "X's"  were  placed  at  each  matrix  posi- 
tion that  corresponded  to  a similarity  of  spectral  signatures  between 
sites.  Note  that  the  CCT  class  limits  in  Table  7 for  sites  1 and  k 
indicate  overlap  on  every  band  and  hence  are  similar.  Thus,  an  X was 
placed  in  the  cluster  matrix  at  the  (l,4)  and  (^+,l)  positions.  Ob- 
viously, the  diagonal  elements  (e.g.  (l,l),  (2,2),  etc.)  all  contain 
X's  since  each  site  is  always  similar  to  itself.  The  procedure  was 
continued  until  every  permutation  of  pairs  of  sites  had  been  compared. 

83.  The  problem  at  this  point  is  to  determine  whether  there  is 
any  order  in  the  array  of  X's  (Block  12.7,  Plate  2).  The  procedure  in- 
volves the  rearrangement  of  the  sequence  of  site  numbers  so  that  simi- 
larities are  brought  together  into  related  "clusters."  The  process  of 
rearrangement  is  called  "diagonalization , " and  results  in  diagonalizing 
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the  cluster  matrix  (i.e.  rearranging  the  cluster  matrix  into  diagonal 
form).  For  example,  site  1 is  similar  to  site  U,  so  the  matrix  is  re- 
arranged to  put  sites  1 and  k into  juxtaposition.  Sites  10  and  12  are 
also  similar  to  site  1,  so  they  too  are  hrought  into  juxtaposition.  The 
result  is  a cluster  matrix  as  illustrated  in  Figure  37. 


SITE  NUMBERS 

I 4 10  12  2 8 15  18  19  20  21  22  2 3 25  27  28  30 


Figure  37.  ClustLM'  matrix,  13  uetober  1975 
data,  with  sites  :,  •,  , aii  i I.’  arranged 

in  a , a.5t‘ a' 

83.  Considering  the  rearranged  matrix,  it  is  clear  tiiat  sites  1, 
10,  and  12  are  not  similar  to  any  -.’ther  .lit*'.-,  and  tiiey  can  therefore 
constitute  a cluster.  Site  2 is  similar  to  site  28,  and  only  28.  Ke- 
arrangement  to  bring  these  two  into  juxtaposition  results  in  the  matrix 
in  Figure  38.  From  the  second  rearrangement  (Figui’e  38),  it  can  be  seen 
that  sites  23  and  25  ai*e  similar  to  site  8,  and  ttiey  can  then  be  brouglit 
together  (Figui'e  39).  It  is  now  evident  that  sites  I8,  27,  and  30  are 
similar  to  site  15.  Bringing  these  together  yields  the  third  arrange- 
ment (Figure  hO) . Examination  of  this  new  arrangement  reveals  an  in- 
complete cluster  formed  by  sites  15,  I8,  27,  and  30.  It  also  discloses 
that  site  22  "overlaps"  the  cluster,  being  similar  to  three  of  the  four 
elements  in  the  cluster.  Moving  site  22  into  juxtaposition  with  the 
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SfTE  NUMBERS 

1 4 10  12  2 28  0 15  18  !9  20  21  22  23  25  27  30 


Figure  33.  Cluster  matrix,  13  October  1975 
data,  with  sites  2 and  28  brought  into  con- 
tiguity 


Figure  39.  Cluster  matrix,  13  October  1975 
data,  with  sites  8,  23,  and  25  brought  into 
contiguity 


SITE  NUMBERS 

I 4 10  12  2 28  8 2 3 2 5 1 5 1 8 27  30  19  20  21  22 


Figure  1+0.  Cluster  matrix,  13  October  1975 
data,  with  sites  15,  l8,  27,  and  30  brought 
into  contiguity 

cluster  yields  the  fourth  arrangement  (Figure  Ul).  This  process  is 
called  diagonalization  because  any  cluster  that  results  is  arranged  along 
the  diagonal,  as  in  Figure  Ul.  Given  a matrix  of  this  form,  a cluster 
can  be  defined  as  a collection  of  elements  (i.e.  occupied  matrix  posi- 
tions) that  can  be  internally  (within  the  cluster)  connected  by  off- 
diagonal  connections  but  that  have  no  connections  with  external  (outside 
the  cluster)  elements  except  by  the  main  diagonal  of  the  cluster  matrix. 

8i+.  The  matrix  in  Figure  !+l  was  interpreted  as  consisting  of  four 
clusters,  three  distinct  and  unique  and  one  somewhat  indefinite  and 
incomplete.  However,  the  fact  that  not  all  matrix  positions  in  the 
cluster  are  filled  does  not  lessen  the  utility.  The  point  is  that  all 
of  the  elements  of  the  incomplete  cluster  are  related  to  each  other, 
even  if  only  indirectly,  while  none  of  its  elements  are  related  to  ele- 
ments in  any  of  the  other  clusters.  Fui’ther,  no  arrangement  of  its  ele- 
ments will  yield  a configuration  of  two  or  more  discrete  nonrelated 
clusters . 


85*  Of  course,  the  site  numbers  can  be  arranged  in  a different 
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SITE  NUMBERS 

1 4 to  12  2 28  8 23  25  15  10  27  30  22  t 9 20  21 
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Figure  i+1.  Cluster  matrix,  13  October  1975 
data;  one  of  several  possible  final  arrange- 
ments showing  four  discrete  clusters 


sequence  (Figure  1*2).  However,  note  that  exactly  the  same  four  clusters, 
one  of  which  is  incomplete,  emerges.  Thus,  the  precise  rationale  by 
which  a site  number  sequence  is  generated  is  unimportant. 

86.  In  the  example  just  described,  fortune  has  smiled:  the 
diagonal! zat ion  process  has  resulted  in  the  formation  of  a cluster 
matrix  containing  discrete  clusters.  However,  as  will  later  be  seen, 
this  is  not  always  the  case.  The  procediire  (Block  12.8,  Plate  2)  for 
dealing  with  cases  in  which  discrete  clusters  cannot  be  formed  is  dis- 
cussed in  paragraphs  110-112. 

8T.  The  initial  terrain  classification  was  to  consist  of  only 
three  categories;  water,  forest,  and  nonforest  (or  open)  (see  para- 
graph 72) • Thus,  each  of  the  sites  in  the  cluster  matrix  (Figure  Ul) 
was  examined  to  determine  whether  the  sites  in  each  cluster  represented 
a single  terrain  type  (Block  12.9,  Plate  2).  By  reference  to  Table  1, 
Table  8 can  be  readily  compiled.  The  comparison  of  Table  8 with  Fig- 
ure !*1  reveals  that  the  four  forest  sites  have  been  neatly  placed  into 
one  of  the  clusters  and  the  open  sites  are  distributed  among  the 
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SITE  NUMBERS 

10  4 1 12  23  25  8 20  21  22  19  30  18  15  27  2 28 


Figure  h2.  Cluster  matrix,  13  October  1975  data,  with 
alternative  arrangement  of  site  numbers  yielding  same 
four  clusters  as  in  Figure  bl 

remaining  three  clusters  (Table  9)-  Each  cluster  then  consists  of  sites 
representing  a sr’ngle  terrain  type  (Block  12. 9 > Plate  2).  If  more  than 
one  terrain  type  had  been  present  in  any  one  of  the  clusters,  the  pro- 
cedure described  in  Blocks  12.20,  12.21,  and  12.22  in  Plate  2 would 
have  been  applied.  No  case  in  this  study  required  this  alternative. 

One  of  the  nonforest  terrain  types,  however,  is  designated  by  more  than 
one  cluster  (Block  12.10,  Plate  2).  The  reason  for  this  can  be  visua- 
lized more  readily  if  the  COT  class  limits  bounding  each  cluster  in  each 
wavelength  band  (Block  12.11,  Plate  2)  are  reduced  to  grapliic  form, 
(Figure  i*3).  The  CCT  values  for  the  limits  bounding  each  cluster  are 
derived  by  determining  the  largest  and  smallest  value  in  each  wavelength 
band  in  each  cluster.  For  example.  Table  10  contains  only  those  sites 
used  to  develop  the  cluster  matrix  in  Figure  h2,  with  the  sites  se- 
quenced in  the  same  order  as  in  Figure  4l.  In  band  4,  since  the 
smallest  value  is  l4  and  the  largest  is  l8,  the  class  limits  of  cluster 
I are  l4-l8.  This  is  plotted  as  a horizontal  bar  in  Figure  43,  as  are 
all  of  the  other  cluster  classes. 
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CCT  UNITS 


Figure  43.  Graphic  representation  of  CCT  classes  derived  from  cluster 
matrix  (Figure  42),  13  October  1975  data,  with  screening  window  and  vari- 
ance window  of  j^l.5  CCT  units 

88.  The  nonforest  clusters  (i.e.  II,  III,  and  IV)  were  grouped 
separately  because  each  failed  the  similarity  tests  in  at  least  one 
spectral  band.  Clearly,  it  would  be  desirable  to  find  a method  of 
arranging  the  sites  such  that  there  is  one  cluster  per  terrain  type.  At 
least  two  ways  exist  in  which  this  might  be  accomplished:  the  size  of 
the  variance  window  might  be  changed  (Block  12.23,  Plate  2),  or  the 
size  of  the  screening  window  might  be  changed  (Block  12.24,  Plate  2),  or 
both  might  be  changed. 

89.  If  the  screening  window  is  reduced,  the  number  of  selected 
sites  will  also  be  lessened  (Table  6).  On  the  other  hand,  if  it  is  in- 
creased, sites  will  be  added  with  the  possibility  that  one  of  them  will 
exhibit  a CCT  value  in  the  critical  wavelength  band  that  will  fill  the 
gap  between  clusters  (Figure  42).  In  the  example,  the  analyst  would  look 
for  a valn.e  that  would  fill  the  critical  gap  between  the  class  limits  of 
clusters  II,  III,  and  IV  in  bands  5 and  6.  However,  note  that  the  class 
limits  of  clusters  II  and  III  overlap  only  band  6.  Thus,  ideally,  the 
sites  that  are  added  should  also  incorporate  values  to  fill  these  gaps. 

90.  From  Table  6,  it  is  apparent  that  widening  the  screening 
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window  to  +1.6  CCT  units  will  add  two  sites,  9 and  13,  to  the  selected 
sites.  The  difficulty  is  that  a reference  to  Table  indicates  that  both 
sites  9 and  13  represent  forest  and  thus  will  serve  only  to  confuse  the 
issue,  i.e.  the  need  was  for  more  sites  in  the  nonforest  areas.  Hence, 
enlarging  the  screening  window  (Block  12.2k,  Plate  2)  will  not  help. 

91.  The  second  possibility  is  to  change  the  variance  window 
(Block  12.23,  Plate  2).  If  the  variance  window  of  those  sites  used  in 
the  cluster  matrix  (Figure  k2)  is  reduced,  the  minimum  effect  will  be  to 
make  the  width  of  the  cluster  classes  smaller,  which  is  opposite  to  the 
effect  desired.  Table  11  illustrates  the  effect  on  the  classes  within 
each  wavelength  band  if  the  variance  window  is  enlarged  to  ^1.7  CCT 
units.  After  a similarity  analysis  is  completed,  a cluster  matrix  is 
assembled,  the  clusters  are  developed  by  the  diagonalization  procedure, 
and  the  cluster  matrix  in  Figure  4k  emerges.  Although  cluster  C is  a 
greatly  extended  cluster,  it  nevertheless  meets  all  of  the  criteria  for 
a cluster  and  also  places  all  nonforest  sample  sites  into  a single 
cluster.  Thus,  the  objective  is  satisfied.  Table  12  is  the  result  of 
arranging  the  class  limits  data  with  the  associated  clusters. 

Site  numbers 

24  10  4 1 12  23  25  8 20  21  22  19  30  18  15  27  2 28 


Fig.ure  44.  Cluster  matrix,  13  October  1975 
data,  with  clusters  derived  by  using,  vari- 
ance window  of  +1.7  CCT  unit.s 
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92.  The  water  site  (site  2h,  Figure  33)  did  not  survive  the  nor- 
malcy analysis  in  the  October  1975  data.  Because  water  so  strongly  ab- 
sorbs band-7  radiation,  the  CCT  value  of  site  2k  in  band  7 should  be  so 
low  that  it  will  not  be  similar  to  any  other  site.  The  data  in  Table  5 
indicate  that  the  mean  CCT  value  of  site  2k  in  band  7 is  2.2,  well  below 
any  of  the  values  exhibited  by  any  other  site,  and  thus  this  expectation 
is  realized.  Since  including  it  among  the  selected  sites  cannot  disturb 
the  normalization  process  (even  though  it  does  not  meet  the  normalcy 
criterion),  site  2k  was  arbitrarily  included  in  the  diagonalized  cluster 
matrix  (Figure  kk)  as  an  independent  one-element  cluster.  The  water 
cluster  could  have  been  added  to  Figure  Ul  in  the  same  way,  but  it  was 
deferred  to  this  point  in  the  procedure  so  that  the  discussion  in  para- 
graph 87  would  be  less  encumbered. 

93.  The  next  procedural  step  was  to  determine  whether  the  terrain 
classifications  were  to  be  "inclusive"  or  "exclusive."  Inclusive  classi- 
fications are  those  in  which  the  classes  include  all  possible  contin- 
gencies; every  patch  of  the  ground  must  be  placed  in  one  or  another  of 
the  categories.  Therefore,  the  analyst  at  this  point  must  decide 
(Block  12.12,  Plate  2)  which  type  is  required.  In  the  present  instance, 
the  system  is  inclusive:  each  patch  of  ground  must  be  either  water, 
forest,  or  nonforest.  Exclusive  classifications  are  those  in  which  only 
certain  aspects  of  a landscape  are  included  in  the  classification,  and 
thus  a certain  proportion  of  the  ground  will  i-emain  outside  the  class 
criteria;  some  patches  will  be  unclassified.  In  such  a case,  the  clas- 
sification actually  consists  of  the  specified  classes  plus  an  additional 
class  in  which  all  things  not  Included  in  the  criteria  defining  the 
specified  classes  are  placed  (Block  12.17,  Plate  2). 

9^*.  The  example  cluster  class  limits  presented  in  Table  12  are 
an  exclusive  set.  However,  since  the  cla.jsification  system  under 
discussion  is  an  inclusive  system  (the  three  clusters  represent  ttiree 
terrain  types  that,  by  definition,  must  cover  100  percent  of  '.lie  ai-ea, 
i.e.,  each  point  on  the  ground  must  lie  either  water,  forest,  i:)r  rriii- 
forest),  the  question  at  this  point  is,  how  are  inclusive  cluster  class 
sets  derived  from  the  exclusive  sets  in  Table  127 

ir) 


95.  The  procedure  (Block  12.13,  Plate  2)  adopted  in  this  study 
determines  by  a "cluster  separation  logic"  the  combinations  of  bands 
that  affect  the  separation  of  clusters  and  then  derives  by  a set  of 
"cluster  group  stretching  rules"  the  inclusive  cluster  class  sets  that 
represent  and  preserve  the  integrity  of  each  cluster  of  the  cluster 
matrix.  A consequence  of  this  proced\ire  is  that  often  more  inclusive 
cluster  class  sets  can  be  determined  than  there  are  original  clusters. 
Furthermore,  since  there  is  not  enough  information  available  to  reliably 
select  the  optimum  set  by  objective  procedures,  a subjective  decision- 
making process  had  to  be  used. 

96.  From  the  graphic  form  of  the  Table  12  data  presented  in 
Figure  U5,  the  cluster  separation  logic  used  to  determine  combinations 
of  bands  affecting  cluster  separations  for  the  October  1975  data  and 
the  symbolic  representation  for  that  logic  is  as  follows : 

a^.  Separation  1:  Clusters  A,  B,  and  C can  be  separated  by 
using  band  6.  The  separation  logic  for  this  cluster 
group  is  represented  by  the  symbol  (Ag  Bg  C). 


CCT  UNITS 

0 5 10  IS  20  2S  30  35  40 


Figure  U5.  Ranges  of  values  comprising  clusters,  13  October  1975 
data,  with  variance  window  of  +1.7  CCT  units 


TO 


b.  Separation  2:  Cluster  A can  be  separated  from  clusters 
B and  C by  using  band  7*  Clusters  B and  C can  be 
separated  from  each  other  by  using  band  5.  The  separa- 
tion logic  for  these  two  cluster  groups  is  represented 
by  the  symbol  (A)  | (B^C). 

c_.  Separation  3:  Cluster  A can  be  separated  from  clusters 
B and  C by  using  band  J.  Clusters  B and  C can  be 
separated  from  each  other  by  using  band  6.  The  repre- 
sentation for  this  separation  logic  is  (A)  I (bI-C). 

( o 

d.  Separation  U:  Cluster  A can  be  separated  from  clusters 
B and  C by  using  band  6,  Clusters  B and  C can  be  sepa- 
rated from  each  other  by  using  band  5-  The  representa- 
tion for  this  separation  logic  is  (a)  g (B^C). 

97.  Once  the  cluster  separation  logics  had  been  determined,  the 
class  limits  of  the  clusters  were  expanded  (or  stretched)  by  means  of  a 
set  of  cluster  group  stretching  rules  (Figure  U6)  to  include  all  possi- 
ble values  in  all  wavelength  bands.  The  rules  are  divided  into  three 
sets:  an  omnigroup,  intragroup,  and  intergroup  set  of  rules.  The  omni- 
group rules  apply  to  all  clusters  regardless  of  their  cluster  groupings; 
the  intragroup  rules,  only  to  bands  affecting  separation  within  a 
cluster  group;  and  the  intergroup  r\iles,  only  to  bands  affecting  a sep- 
aration between  cluster  groups. 

98.  Table  13  presents  the  consequences  of  the  cluster  group 
stretching  rules  as  they  were  successively  applied  to  the  cluster  sep- 
aration logic  (a)  I (B^C)  (i.e.  Separation  2,  paragraph  96).  Table  1^4 
illustrates  the  resiilts  of  the  group  stretching  rules  on  all  the  separa- 
tion logics.  At  this  point  in  the  process,  it  was  still  undetermined 
which  of  the  class  sets  (Table  1J4)  would  produce  the  best  map  of  the 
area  of  interest.  In  this  context,  the  best  map  would  be  that  one  in 
which  the  maximum  number  of  pixels  was  correctly  identified  as  to  ter- 
rain type.  In  the  example  used  up  to  this  point,  it  is  extremely  dif- 
ficult to  estimate  which  will  produce  the  best  product,  since  all  of 
them  are  capable  of  clearly  separating  all  three  clusters.  The  only 
obvious  way  of  resolving  the  problem  was  to  actually  produce  prelimin- 
ary maps  (Blocks  12.1^4  and  12.15,  Plate  2)  of  a small  region  and  check 
the  results  against  the  actual  terrain  and  against  an  air  photo  that  was 
reliably  interpreted.  Figures  ^47-50  are  the  maps  resulting  from  using 
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OmnlRroup  Rules 


Rule  1 
Rule  2 

Rule  1 
Rule  2 
Rule  3 

Rule  4 


Rule  1 
Rule  2 
Rule  3 

Rule  4 


The  lower  class  limits  in  bands  not  affecting  cluster  or 
cluster  group  separations  are  extended  to  0. 

The  upper  class  limits  in  bands  not  affecting  cluster  or 
cluster  group  separations  are  extended  to  63. 


Intragroup  Rules 

The  lower  limiting  CCT  value  in  the  band  affecting  separation 
is  extended  to  0. 

The  upper  limiting  CCT  value  in  the  band  affecting  separation 
is  extended  to  63. 

In  the  band  affecting  separation,  the  upper  and  lower  class 
limits  of  a cluster  are  extended  until  the  gaps  between 
adjacent  clusters  are  closed.  If  the  gap  between  an  adjacent 
cluster  is  an  even  number  of  CCT  value  units,  divide  the 
values  equally.  If  the  gap  is  an  odd  number  of  CCT  value 
units,  assign  the  extra  value  to  the  lower  of  the  two  clusters. 
If  the  cluster  group  consists  of  only  one  cluster,  extend  its 
lower  and  upper  class  limits  to  0 and  63,  respectively,  in 
every  band  not  affecting  an  intergroup  separation. 


Intergroup  Rules 

The  lowest  CCT  value  in  the  band  affecting  separation  is 
extended  to  0. 

The  highest  CCT  value  in  the  band  affecting  separation  is 
extended  to  63. 

In  the  band  affecting  separation,  the  lowest  and  highest  CCT 
class  limits  of  a cluster  group  are  assigned  to  each  cluster 
of  that  group.  For  example,  clusters  X,  Y,  and  Z are  in  the 
same  cluster  group  and  have  class  limits  of  21-25,  23-29,  and 
19-24  in  band  6,  respectively;  then  the  band  6 class  limits 
assigned  to  X,  Y,  and  Z are  19-29. 

In  the  band  affecting  separation,  the  upper  and  lower  limits 
of  a cluster  group  are  extended  until  the  gaps  between 
adjacent  cluster  groups  are  closed.  If  the  gap  between  an 
adjacent  cluster  group  is  an  even  number  of  CCT  value  units, 
divide  the  number  equally.  If  the  gap  is  an  odd  number  of 
CCT  value  units,  assign  the  extra  value  to  the  lower  of  the 
two  cluster  groups. 


Fif^ure  46.  Cluster  group  stretching  rules 
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Fif.ure  ^47.  Study  area  aa  interpreted  by  cluater  clasa  set  i ^band  6, 

13  October  1975  data) 
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Study  area  as  interpreted  by  cluster  class  set  2 
(bands  5 and  7,  13  October  1975  data) 


7't 


Figure  U8. 
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Figure  k9.  Study  area  as  interpreted  by  cluster  class  set  3 (bands  6 

and  7,  13  October  1975  data) 
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Fi(-;ure  50.  Study  area  as  interpreted  by  cluster  class  set  U (bands  5 ' 

and  6,  13  October  1975  data)  j 
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the  "stretched"  CCT  class  sets  given  in  Table  lU, 

99-  A comparison  between  the  maps  generated  by  the  four  cluster 
class  sets  and  an  air  photo  of  the  same  region  (Figure  lU)  indicates  that 
cluster  class  set  2 (Figure  k8} , which  employs  bands-5  and  -7  radiance 
data,  provides  the  best  interpretation  of  the  terrain  conditions 
(Block  12. l6,  Plate  2).  The  maps  produced  by  cluster  class  sets  1 and 
3 (Figures  h'J  and  U9)  depict  numerous  patches  of  nonforest  (cluster  C) 
within  the  large  forest  areas.  Also,  set  1 does  not  detect  the  presence 
of  the  river  in  the  upper-left  part  of  the  study  area.  The  map  produced 
by  cluster  class  set  U (Figure  50 ) also  does  not  detect  some  of  the  promi- 
nent water  bodies  (the  river  and  a small  lake  in  the  west-central  part 
of  the  study  area).  On  the  other  hand,  cluster  class  set  i*  did  give  an 
excellent  interpretation  of  the  distribution  of  forest  and  nonforest  in 
the  study  area. 

100.  The  subjective  impressions  described  above  were  then  checked 
by  comparing  the  area  of  each  class,  as  determined  by  automatic  classi- 
fication, with  the  areas  as  carefully  mapped  by  air-photo  interpretation 
by  a skilled  interpreter  with  excellent  ground  truth  data  (Figure  lU). 

The  areas  in  Figure  lU  were  carefully  measured  using  a dot-count  tech- 
nique. Table  15  gives  the  resultant  comparisons.  The  results  are  some- 
what surprising.  Note  that  the  errors  from  set  3 are  somewhat  less  than 
those  obtained  from  set  2,  despite  the  fact  that  the  subjective  impres- 
sion is  that  the  map  produced  by  set  2 is  superior  to  that  produced  by 
set  3.  The  question  is,  why  does  not  the  subjective  impression  match 
the  objective  analysis?  The  apparent  reason  is  subtle.  The  areas 
presented  in  Table  15  are  meaningful  only  as  totals;  there  is  no  neces- 
sary relation  between  the  total  area  covered  by  a category  and  the  lo- 
cation of  individual  pixels  of  that  category.  Thus,  at  least  in  princi- 
ple, every  pixel  could  be  classified  wrongly,  and  yet  the  totaled  areas 
could  agree  perfectly.  The  obvious  conclusion  is  that  total  areas  are 
not  meaningful  as  a basis  for  a judgment  as  to  which  of  several  classi- 
fication schemes  is  superior.  The  only  valid  basis  is  to  make  the 
comparison  on  a pixel-by -pixel  basis.  When  this  is  done  (Figures  lU  and 
^T-50),  it  will  be  noted  that  cluster  class  sets  1 and  3 have  many 
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obviously  misclassif ied  pixels  scattered  throughout  both  the  forest  and 
nonforest  areas.  The  remarkable  agreement  between  total  areas  is 
fortuitous.  Cluster  class  set  h has  far  fewer  misclassified  pixels  in 
the  forest  and  nonforest  areas  and  is  deficient  only  in  the  fact  that  it 
does  not  appear  to  map  water  surfaces  as  effectively  as  set  2.  However, 
on  balance,  cluster  class  set  2 was  chosen  as  the  best  of  the  classifi- 
cation schemes  for  the  October  1975  data. 

101.  In  the  event  that  the  classification  system  is  an  exclusive 
type  (paragrapn  93),  the  procedure  is  to  go  directly  from  Block  12.17 
to  Block  13  (Plate  l).  In  this  case,  of  course,  the  cluster  class 
sets  will  include  not  only  class  limits  for  the  desired  terrain  types 
but  also  an  "unclassified"  category,  which  will  accumulate  all  pixels 
not  meeting  the  criteria  of  the  cluster  class  sets. 

102.  Block  13.  Having  made  all  of  the  decisions  required  to 
classify  every  pixel  in  the  scene  or  area  of  concern,  the  procedure  was 

to  do  exactly  that,  using  the  original  geometrically  uncorrected  CCT's.  ! 

The  program  that  performed  this  operation  accepts  the  limiting  CCT  class 
values  for  each  cluster  (which  represents  a terrain  type)  in  each  wave- 
length band.  The  four  CCT  values  defining  each  pixel  are  then  compared 

with  the  CCT  class  specifications  (Table  lU),  and  that  pixel  is  then  | 

coded  with  a number  that  specifies  the  terrain  type  that  correlates  with  j 

the  set  of  CCT  classes  into  which  it  fits.  The  result  is  a new  CCT  con-  i 

taining  only  the  category  codes  instead  of  the  original  sets  of  CCT  j 

values . 

103.  Block  lU.  The  next  procedural  step  was  to  bring  the  new  j 

CCT  containing  only  the  category  codes  into  geometric  accordance  with  ( 

t 

the  air  photos.  The  Landsat  pixel  array  could  have  been  brought  into  | 

geometric  accordance  prior  to  classification.  There  are,  however,  j 

compelling  reasons  for  not  doing  so.  The  most  important  one  revolves  j 

around  the  matter  of  pixel  shape.  The  Landsat  pixel  is  rectangular,  j 

approximately  58.1  by  79-2  m;  the  pixel  produced  by  the  Optronics  film 
reader /writer  (see  Appendix  A)  is  a perfect  square.  The  discrepancy  can 
be  approximately  corrected  by  compiling  a new  pixel  array  in  which 
every  3rd  and  20th  scan  line  of  the  original  CCT  is  duplicated.  This 
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procedure  increases  tne  number  of  scan  lines  in  a Landsat  scene  from 
23UO  to  3237 • Thus,  there  are  approximately  1.36  times  as  many  pixels 
in  the  reformatted  CCT  as  in  the  original  Landsat  CCT.  The  new  tape 
can  now  be  used  on  the  film  reader/writer  to  produce  an  image  that  is 
approximately  the  same  scale  both  parallel  to  and  at  right  angles  to 
the  orbital  path;  whereas  had  the  original  tape  been  used,  the  film 
reader /writer  would  have  produced  an  image  in  which  the  scales  differed 
by  a factor  of  approximately  1:1.36. 

10i+.  The  important  aspect  of  this,  with  respect  to  the  classi- 
fication process,  is  that  all  of  the  available  data  are  in  the  original 
CCT;  reformatting  adds  pixels  but  not  information.  Thus,  if  all  manipu- 
lations required  in  the  classification  procedure  are  done  prior  to  re- 
formatting, substantial  savings  in  computer  time  are  effected,  simply 
because  fewer  pixels  need  be  handled. 

105.  The  second  major  reason  for  deferring  shape  correction  stems 
from  the  image  skew  produced  by  the  rotating  earth.  Approximately 
2h  sec  are  required  for  Landsat  to  move  I85  km  (the  length  of  a Landsat 
scene)  along  its  orbit,  and  during  this  time,  the  rotation  of  the  earth 
carries  the  surface  eastward  about  11.49  km  at  the  latitude  of  Vicksburg, 
Mississippi.  This  means  that  the  surface  moves  eastward  with  respect  to 
ohe  satellite  orbital  path  one  pixel  width  (57  m)  in  about  0.124  sec. 
Therefore,  in  the  time  that  it  takes  the  surface  to  move  eastward  the 
width  of  one  pixel,  the  satellite  moves  southward  along  the  orbit  about 
11.6  scan  lines.  This  systematic  error  can  be  approximately  corrected 
(see  Reference  8 for  a more  accurate  correction)  by  offsetting  each 
successive  group  of  12  scan  lines  to  the  westward  one  pixel  widtli.  For 
an  area  the  size  of  a Lar.isat  scene,  a total  of  260  "false  pixels"  are 
inserted  on  the  tape,  divided  appropriately  between  the  ends  of  each 
scan  line,  to  bring  the  digital  array  back  into  a rectangular  array. 

The  result  of  this  process  is  only  barely  perceptible  at  a scale  of 
1:1,000,000  on  the  resulting  film  because  of  the  small  pixel  size  used 
to  expose  the  film,  but  it  is  nevertheless  effective  for  approximately 
correcting  for  skew.  The  effect,  however,  is  to  add  841,620  pixels  to 
the  data  array.  If  these  pixels  are  included  in  the  data  array  prior 
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to  classification,  special  arrangements  must  be  made  so  that  they  are 
not  included  among  those  pixels  to  be  classified.  Regardless  of  how 
this  is  done,  the  process  will  require  computer  time.  Thus,  conser- 
vation of  computer  time  is  achieved  by  deferring  geometric  correction 
^intil  after  the  classification  procedure  is  completed. 

106.  Block  15.  The  final  step  is  to  write  a map  showing  the 
distribution  of  terrain  types,  as  determined  by  the  procedure  described 
above.  The  point  of  departure  is  the  reformatted  CCT  in  which  each 
pixel  is  identified  by  a "terrain  type"  code  (see  paragraph  102)  and  in 
which  each  pixel  is  located  in  such  a way  that  it  will  be  located  in  the 
proper  position  when  fitted  over  a base  map. 

107.  The  light-emitting  diode  (LED)  in  the  Optronics  film  reader/ 
writer  (see  Appendix  A)  is  capable  of  producing  256  light-intensity 
levels,  which  results  in  as  many  "gray  shades"  on  the  exposed  film. 

Since  the  human  eye  has  difficulty  distinguishing  more  than  about  five 
to  seven  shades  of  gray,  it  is  usually  wise  to  restrict  the  number  of 
categories  on  any  one  map  to  no  more  than  seven  different  categories. 

In  the  present  instance,  since  only  three  categories  were  required,  this 
constraint  offered  no  difficulties.  The  gray  shades  that  were  selected 
for  best  portrayal  of  the  three  categories  were  as  follows : 


Mapping  Cluster 

Terrain  Type 

LED  Intensity 
Level 

1 

Gray  Shade  1 

A 

Water 

255 

1 

Black 

B 

Forest 

70 

Gray 

C 

Nonforest 

0 

White 

These  values  were  used  to  produce  a map  of  the  study  area  (Figure  51 ), 
and  also  a map  of  the  entire  October  1975  Landsat  scene  (Figure  52), 
each  classifying  the  areas  with  the  three  terrain  types. 

Supplemental  tests  of  procedure 

108.  Since  only  one  test  of  a procedure  constitutes  a somewhat 
uncertain  basis  for  evaluating  its  utility,  two  additional  Landsat  scenes 
(Table  2)  covering  essentially  the  same  region  were  processed.  The 
differences  in  season  (note  that  one  scone  was  obtained  in  the  winter 


80 


i 


4 


(21  February  1975),  one  in  summer  (11  July  197^),  and  another  in  the 
autumn  (13  October  1975))  were  expected  to  produce  significant  changes 
in  terrain  reflectance  properties  and,  therefore,  in  spectral  signatures. 
The  sun  position  also  chaiages  significantly  and  thus  causes  changes  in 
spectral  signatures  because  of  differences  in  reflectance  geometry. 
Finally,  normal  agricultural  practices  could  be  expected  to  change 
field  characteristics  in  the  cultivated  areas.  These  factors  all 
conspire  to  make  certain  that  the  cluster  class  sets  derived  for  the 
October  1975  scene  would  not  be  suitable  for  the  other  two  scenes.  In 
effect,  the  analytical  procedure  must  start  from  the  beginning  for  each 
scene . 

109.  However,  the  same  test  area  (Figures  U and  5),  as  well  as 
the  same  sample  sites  (Figure  33)  was  used  in  all  cases.  Further,  a 
modest  amount  of  ground  truth  data  had  been  collected  (see  Part  III) 
prior  to  the  analysis  of  the  imagery  data,  and  thus  the  terrain  type  at 
each  sample  site  was  known  (Table  4).  With  this  somewhat  tenuous  ground 
truth  data,  the  July  1974  and  February  1975  CCT's  were  subjected  to 
exactly  the  same  procedure  as  the  October  1975  COT. 

110.  Analysis  of  11  July  1974  data.  Table  16  presents  the  result 
of  the  normalcy  analysis  on  the  3-  by  3-pixel  arrays  around  the  30  sam- 
ple sites.  Table  17  shows  the  effects  of  changes  in  the  size  of  the 
screening  window  on  the  number  of  sample  sites.  For  instance,  only  10 
sample  sites  survive  with  a screening  window  of  +1.5  OCT  units,  and 
only  15  survive  with  a screening  window  of  +^1.8  COT  units.  As  a general 
principle,  the  degree  of  representativeness  can  be  expected  to  be  di- 
rectly related  to  the  number  of  sample  sites.  The  initial  cluster 
matrix  was  constructed  with  those  sample  sites  derived  from  a screening 
window  of  j^l.8  COT  units  and  a variance  window  (Table  I8)  of  +1.5  OCT 
units.  After  diagonal izat ion , the  cluster  matrix  is  arranged,  as  in 
Figure  53.  At  this  point,  the  procedure  is  at  Block  12.8,  Plate  2,  and 
it  is  evident  that  clean  clusters  cannot  be  formed.  There  is  a possi- 
bility that  a smaller  variance  window  will  eliminate  the  cluster  overlap 
(Block  12.18,  Plate  2). 
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Figure  53.  Cluster  matrix,  11  July  197^  data, 
with  screening  window  of  Hy..8  CCT  units  and 
variance  window  of  ^1.5  CCT  units 

111.  Reducing  the  variance  window  to  j^l.3  CCT  units  (the  smallest 
acceptable  window)  results  in  the  cluster  matrix  in  Figure  5^.  It  is 
clear  that  clean  clusters  cannot  be  formed.  One  alternative  remains, 
namely  that  of  reducing  the  screening  window  (Block  12.19,  Plate  2),  in 
the  hope  that  clean  clusters  will  form. 

112.  A screening  window  of  ^1.5  CCT  units  was  then  selected,  pri- 
marily because  it  was  successful  for  the  October  1975  data.  Figure  55 
shows  the  cluster  matrix  that  results  from  using  the  variance  window  of 
j^l.5.  Since  the  clusters  have  not  yet  cleanly  formed,  the  obvious  thing 
to  do  would  be  to  reduce  the  variance  window.  Tlie  use  of  a smaller 
variance  window  (^1.3  CCT  units)  results  in  the  CCT  classes  given  in 
Table  19  and  the  cluster  matrix  in  P’igure  56.  Figure  57  is  a graphic 
presentation  of  the  resulting  cluster  class  limits. 

113.  An  examination  of  Figure  57  indicates  that,  just  as  with 
the  October  1975  data,  separation  of  water  is  very  easy  with  either 
band  6 or  band  7.  Although  previous  experience  had  indicated  that 
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Figure  5i+.  Cluster  matrix,  11  July  197^  data, 
with  screening  window  of  j^l.8  CCT  units  and 
variance  window  of  +1.3  CCT  units- 


Fif’;ure  55-  Cluster  matrix,  11  July  197^  data, 
with  screening  window  and  variance  window 
of  +1.5  CCT  units 
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Figure  56.  Cluster  matrix,  11  July  197^  data, 
with  screening  window  of  j^l.5  CCT  units  and 
variance  window  of  +1.3  CCT  units 
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Figure  57-  Graphic  display  of  cluster  class  limits,  11  July  197^  data, 
with  screening  window  of  • 5 CCT  units  and  variance  window  of  ^1.3 

CCT  units 
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band  7 would  yield  superior  results,  it  was  decided  to  test  both  bands. 

Clusters  B and  C (forest  and  nonforest)  can  obviously  only  be  cleanly 
separated  by  the  use  of  band-5  data.  Table  20  presents  the  resulting 
stretched  cluster  class  sets  as  eventually  used.  Figure  58  shows  the 
final  map  product  of  the  use  of  cluster  class  set  1. 

11^.  Analysis  of  21  February  1975  data.  The  same  sample  sites 
(Figure  33)  were  used  for  the  analysis  of  the  February  1975  as  for  the 

July  197^  and  October  1975  data  sets.  Table  21  gives  the  mean  and  j 

standard  deviations  of  the  signatures  taken  from  the  sample  3-  by  3-pixel 
arrays,  and  Table  17  lists  the  sites  that  survived  normalcy  evaluation 
at  various  screening  windows.  Table  22  shows  the  effect  of  changes  in 
the  variance  window  on  the  CCT  class  limits  of  the  17  sites  that  sur- 
vived normalcy  analysis  with  a screening  window  of  +1.8  CCT  units. 

115.  Figure  59  is  the  cluster  matrix  that  results  from  the  use 
of  a variance  window  of  +1.5  CCT  units,  and  Figure  60  is  the  matrix  in 
diagonal  form.  It  is  evident  that  only  two  clusters  can  be  fox’med,  one 
consisting  of  only  one  site  (site  29).  Reducing  the  variance  window  to 
j^l.3  CCT  units  yields  a cluster  matrix,  as  shown  in  Figure  6l.  It  is 
clear  that  site  13  is  the  anomalous  site;  its  removal  results  in  a 
clean  separation  between  forest  and  nonforest,  as  indicated  by  the  site 
characteristics  (Table  U)  along  the  right  margin  of  the  figure.  If  it 
co\ild  be  removed  from  the  set  of  site  numbers,  all  would  be  well. 

116.  The  procedure  for  removing  site  13  from  contention  is  to 

reduce  the  screening  window.  Table  17  indicates  that  site  13  will  be 
eliminated  if  the  screening  window  is  reduced  to  CCT  units.  Fig- 

ure 62  shows  the  cluster  matrix  resulting  from  the  use  of  a screening 
window  and  a variance  window  of  ^1.5  CCT  units.  Though  the  forest  sites 

form  two  discrete  clusters  and  are  somewhat  underrepresented  (2  out  of  | 

12  sites),  at  least  the  nonforest  sites  now  form  a discrete  cluster.  j 

Figure  63  is  a graphic  representation  cf  the  resulting  cluster  class  | 

limits  (Table  23).  The  water  site,  site  2h , has  been  arbitrarily  added,  i 

even  though  it  did  not  meet  the  normalcy  criterion.  | 

117.  As  noted  in  Figure  63,  band  7 can  be  used  to  separate  j 

clusters  A and  from  clusters  and  C;  band  5 offers  an  opportunity  | 
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Figure  59-  Cluster  matrix,  21  February  1975 
data,  with  screening  window  of  +1.8  CCT  units 
and  variance  window  of  +1-5  CCT  units 
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Figure  60.  Diagonalized  cluster  matrix, 
21  February  1975  data,  with  screening 
window  of  jy..8  CCT  units  and  variance 
window  of  +1.5  CCT  units 
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Figure  63.  Graphic  presentation  of  cluster  class  limits,  21  February 
1975  data,  with  screening  window  and  variance  window  of  jy-.5  CCT  units 

of  sharply  separating  cluster  A from  cluster  and  cluster  from 
cluster  C;  and  band  6 can  accomplish  the  same  separations  as  band  7- 
Thus,  two  possible  cluster  class  sets  exist  that  will  separate  the  four 
clusters. 

118.  With  the  cluster  class  distributions  (Figure  63)  given,  two 
stretched  cluster  class  sets  were  determined  (Table  2h) . A test  indi- 
cated that  cluster  class  set  2 was  superior.  The  map  of  the  study  area 
produced  with  these  cluster  class  limits  (Figure  61+)  shows  the  two 
forest  clusters  B^  and  B^  grouped  together  and  displayed  as  one  shade  of 
gray. 

119.  The  formation  of  the  two  forest  clusters,  B^^  and  B^,  rather 
than  just  one  discrete  cluster  can  be  understood  by  referring  to  the 
cluster  class  limits  in  Table  2k.  It  is  evident  from  the  table  that 
signatures  from  Landsat  pixels  belonging  to  the  forest  cluster  cannot 
be  grouped  with  the  "normal"  forest  cluster  B^,  because  their  band-7  CCT 
values  are  too  low  and  very  near  the  values  of  water.  On  the  other  hand, 
they  cannot  be  grouped  with  water  either,  since  their  band-5  CCT  values 
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are  much  lower  than  those  of  water.  The  final  interpretation  of  this 
situation  was  that  the  pixels  that  fit  the  cluster  class  set  were  a 
mixture  of  water  and  forest  terrain  types  and  thus  formed  naturally  into 
a separate  and  discrete  cluster.  This  is  discussed  further  in  greater 
detail  in  paragraphs  13^-136. 
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PART  III:  EVALUATION  AND  ADJUSTMENTS 


Evaluation 


120.  Any  semiautomated  terrain  interpretation  process  using 

Landsat  digital  data  will  almost  inevitaUly  produce  "patches"  of  mis- 

classification,  because  some  pixels  will  fall  in  locations  such  that 

parts  of  two  or  more  different  terrain  types  are  represented.  The  GOT 

values  of  such  pixels  will  be  a mix  of  the  spectral  signatures  of  all  of 

the  terrain  types  included  in  the  pixel.  Such  "spurious"  signatures  may 

resemble  no  known  terrain  type,  or  they  may,  by  chance,  resemble  a ter- 

1* 

rain  type  other  than  those  actually  encompassed  by  the  pixel. 

121.  There  are  also  other  reasons  for  misclassification.  For 
example,  the  entire  procedure  depends  on  the  assumption  that  different 
terrain  types  will  always  exhibit  different  spectral  characteristics. 
This  assumption  is  not  necessarily  valid.  It  is  entirely  possible  that 
a specific  reflectance  geometry  will  shift  the  spectral  signature  of 
one  terrain  type  to  values  closely  resembling  or  identical  to  those  of 
another  terrain  type.  Furthermore,  two  completely  different  terrain 
types  may  have  the  same  reflectance  properties , in  which  case  it  will  be 
impossible  to  differentiate  them  by  spectral  analysis. 

122.  The  essential  point  is  that  terrain  classification  on  the 
basis  of  spectral  analysis  alone  cannot  be  expected  to  achieve  perfect 
discrimination.  The  goal  is  to  keep  misclassification  to  a minimum. 
Accordingly,  the  final  products  were  sub,;ected  to  close  examination  to 
determine  sources  of  misclassification.  in  the  hope  that  such  informa- 
tion could  be  used  to  improve  the  procedure. 

123.  Unfortunately,  the  mc'^t.  recent  air  photos  of  the  study  area 
that  could  be  found  were  taken  in  mid-July  1973  (Figxire  5).  Thus,  it 
was  impossible  to  compare  the  terrain  type  maps  produced  by  the  semi- 
automated spectral  analysis  procedure  with  actual  ground  conditions  as 
revealed  by  air-photo  interpretation.  However,  it  was  possible  to  visit 
the  study  area  on  the  day  of  the  October  1975  overpass  and  to  record 
the  conditions.  Figures  16-20  illustrate  some  of  ttie  critical 
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conditions,  and  Figure  21  indicates  their  ground  locations.  The  area 
was  visited  again  on  8 May  1976,  and  critical  points  were  photographed 
from  the  air  at  low  altitude.  A set  of  these  photographs  are  included 
as  Fig\ires  22-30,  and  their  locations  indicated  in  Figvire  31.  The 
combination  of  ground  truth  data  and  low-altitude  examination  from  the 
air  provided  detailed  control  on  the  locations  of  critical  boundaries 
and  the  identification  of  critical  areas. 

12l*.  In  addition,  local  people  were  interviewed  in  an  attempt  to 
establish  conditions  at  the  times  of  the  July  197^  and  February  1975 
overpasses.  The  conditions  in  every  patch  could  not  be  recorded,  but 
the  crops  in  and  conditions  of  certain  fields  were  established  with  con- 
siderable reliability.  Figures  12  and  13  illustrate  much  of  this  infor- 
mation. The  major  difficulty  is  that  the  floodplain  was  still  partly 
in\mdated  by  the  water  of  the  great  floods  of  197^-75  (Figure  7)  during 
the  time  of  the  July  197^  and  February  1975  overpasses,  and  this  is,  of 
course,  not  visible  on  the  July  1973  air  photo.  In  an  attempt  to  antic- 
ipate what  might  show  on  the  two  Landsat  images,  two  water  levels,  one 
at  85  ft  to  correspond  to  the  July  197^  overpp""  i,  and  another  at  90  ft 
to  correspond  approximately  to  the  Februar’  overflight,  were  pro- 

jected onto  the  best  available  topographic  jf  the  study  area.  Fig- 
ure 65  shows  the  result.  Note  that  large  areas  of  floodplain  on  the 
western  edge  of  the  study  area  were  subject  to  inundation  at  both  the 
85-  and  90-ft  levels. 

125.  Too  much  credence  should  not  be  placed  in  the  precise  loca- 
tions of  the  limits  of  the  projected  flooded  areas.  There  are  two  major 
reasons.  First,  a very  small  error  in  measuring  elevations  can  move  the 
position  of  a contour  line  a substantial  distance  in  such  flatland  as 
the  Mississippi  floodplain.  The  existing  topographic  maps  are  largely 
the  product  of  photogrammetric  analysis,  and  it  is  extremely  difficult 
to  determine  elevations  accurately  by  that  method  in  regions  of  very  low 
relief.  As  a consequence,  contour  lines  in  the  Mississippi  floodplain 
are  not  uncommonly  misplaced  by  tens  or  even  hundreds  of  metres.  The 
second  reason  is  that  the  elevations  plotted  in  Figure  7 reflect  gage 
readings  on  the  Yazoo  River  at  Satartia.  The  gage  is  very  close  to  the 
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study  area,  as  indicated  in  Figure  1+,  but  the  water  levels  on  the  flood- 
plain  itself  do  not  march  to  exactly  the  same  rhythm  as  the  river.  In- 
stead, they  normally  lag  well  behind  the  river  gages.  Thus,  for  example, 
the  July  7^  Landsat  overpass  occurred  during  a falling  stage  (Figure  7), 
and  the  water  on  the  floodplain  could  be  expected  to  be  at  some  in- 
determinately higher  level  than  the  gage  reading,  simply  because  the 
water  would  not  have  had  time  to  drain  off.  On  the  other  hand,  the 
February  1975  overpass  occurred  on  what  was  essentially  a rising  stage, 
and  thus  the  levels  on  the  floodplain  could  be  expected  to  be  somewhat 
lower  than  the  gage  reading,  because  the  low  areas  would  not  yet  have 
had  time  to  fill.  The  difficulty  is  that  in  neither  case  can  the  exact 
amount  of  lag  be  predicted  with  complete  reliability. 

126.  As  anticipated,  the  most  striking  differences  among  the 
three  final  products  of  the  classification  system,  as  applied  to  the 
test  area,  are  in  the  delineation  of  water  areas  (Figures  51,  58,  and 
6i|).  For  example,  the  July  197^  map  (Figure  58)  shows  very  extensive 
water  areas,  especially  at  the  western  edge  of  the  study  area.  When 
this  is  compared  with  the  map  of  potentially  flooded  areas  (Figure  65), 
it  will  be  seen  that  the  areas  mapped  as  water  coincide  closely  with  the 
areas  of  possible  inundation.  The  absolute  distributions  will  not  match, 
of  course,  for  the  reasons  discussed  in  paragraph  125.  Very  close 
examination  of  the  photomosaic  (Figure  5)  will  reveal  that  the  tiny 
water  areas  scattered  through  the  forest  area  covering  the  southeast 
half  of  the  study  area  and  in  the  cultivated  area  in  the  southeast 
quadrant  represent  small  reservoirs  or  farm  ponds. 

127.  The  map  of  the  February  1975  data  also  shows  a considerable 
amoxint  of  water  along  the  western  margin  of  the  test  area  (Figure  6h), 
but  it  is  much  less  than  that  mapped  from  the  Jxily  197^  data  (Figure  58). 
A comparison  with  the  map  of  potential  flooding  (Figure  65)  shows  that 
the  areas  mapped  as  water  coincide  closely  with  the  areas  below  the  85- 
ft  contour  line.  For  reasons  discussed  in  paragraph  125,  this  is  not 
unreasonable. 

128.  The  map  of  the  study  area  based  on  October  1975  data 
(Figxire  51)  indicates  that  the  water  had  drained  off  the  floodplain. 
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just  as  would  be  anticipated  by  considering  the  Yazoo  River  gage  levels 
(Figure  J).  Note  also  that  the  Yazoo  River,  which  loops  through  the 


extreme  northwest  corner  of  the  test  area,  is  so  narrow  that  it  shows  as 
a discontinuous  line  of  pixels  ( see  Reference  for  a detailed  explana- 
tion of  the  phenomena  responsible);  whereas,  in  the  July  197^  and  February 
1975  scenes  (Figiires  58  and  6h,  respectively),  the  river  is  mapped  as  a 
continuous  line  of  water  pixels. 

129.  It  may  be  concluded  that  the  procedure  provides  a highly 
reliable  method  of  mapping  water  surfaces,  even  when  the  water  is  quite 
turbid  (as  it  is  in  the  Yazoo  River  in  time  of  flood)  and  very  shallow 
(as  it  must  have  been  over  large  areas  of  the  floodplain).  However, 
there  are  some  anomalies.  For  example,  the  large  area  centered  approxi- 
mately at  coordinates  3. 5-7.0  is  classified  as  forest  in  Figures  58  and 
6U,  although  it  ought  to  be  flooded  according  to  the  topographic  map 
(Figure  65). 

130.  The  difficulty  lies  in  the  fact  that  the  two  terrain  types 
are  not  exclusive.  The  area  is  actually  a large  patch  of  floodplain 
forest,  and  ground  examination  revealed  that  the  area  beneath  the  trees 
had  been  inundated  during  the  flood  period.  However,  the  Landsat  sensor 
received  radiation  from  the  forest  canopy  and  not  from  the  water  beneath 
the  trees.  Thus,  the  classification  algorithm  correctly  classified  the 
area  as  forest.  Close  examination  will  reveal  other  similar  examples. 

131.  Some  anomalies  are  not  so  readily  explainable.  For  instance, 
the  February  and  October  1975  maps  (Figures  6k  and  51 » respectively) 
both  indicate  a patch  of  forest  at  about  coordinates  6.0-10.0,  but  the 
July  197^  map  (Figure  58)  shows  none  at  that  location.  The  1973  air 
photo  (Figure  5)  indicates  that  there  was  a patch  of  forest  at  that 
location  in  July  1973.  Why  did  the  algorithm  misclassify  it  as  non- 
forest in  the  July  197^  data?  It  is  not  yet  known. 

132.  However,  the  Landsat  maps  do  reveal  the  removal  of  forest. 

For  example,  the  large  patch  of  forest  centered  on  coordinates  7.6-12.0 
(Figure  51)  is  indicated  on  the  1973  air  photo  (Figure  5)  as  having  a 
long  tongue  extending  to  the  eastward,  but  none  of  the  Landsat  scenes 
show  its  presence.  Ebcamination  on  the  ground  revealed  that  it  had 
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been  cleared  and  converted  to  agricultviral  land  (Figures  lU  and  15b). 

The  same  situation  is  evident  at  the  southeast  corner  of  the  large  patch 
(Figure  15a). 

133.  The  flooded  forest  region  at  coordinates  3.T-7.0  also  con- 
tains an  example  of  misclassification.  In  the  July  197^  map  (Fig- 
ure 58),  there  is  a discontinuous  line  of  nonforest  pixels  along  the 
western  and  southern  edges  of  the  forest  patch.  The  topographic  map 
(Figure  9)  indicates  no  elevation  changes  that  co\ild  bring  this  small 
zone  above  the  level  of  the  floodwaters;  the  area  must  have  been 

1* 

flooded.  The  anomaly  is  probably  the  result  of  spurious  signatures. 

The  pixels  in  this  region  fell  partly  on  forest  and  partly  on  water,  and 
by  chance  the  integration  of  the  radiance  values  produced  a signature 
that  met  the  nonforest  criteria. 

13^.  Why  did  the  same  phenomenon  not  occur  in  the  February  1975 
data  (Figure  6U)?  A reference  to  Figure  63  will  reveal  the  probable 
reason.  Note  that  one  of  the  forest  clusters  (cluster  B^)  exhibits 
CCT  values  in  band  7 so  low  that  they  overlap  the  CCT  values  of  the 
water  cluster  (cluster  C).  Further,  even  in  band  5 (Figure  6U  was 
constructed  with  cluster  class  set  2;  see  Table  2k),  the  separation 
between  water  and  forest  clusters  is  not  large.  Thus,  a spurious 
signature  caused  by  a pixel  falling  across  a boundary  is  far  more  likely 
to  be  classed  as  one  of  the  two  elements  comprising  the  combination  than 
some  third  category. 

Ad.iustments 

135.  The  fact  that  a forest  cluster  looked  so  much  like  a water 
cluster  in  the  February  1975  data  led  to  a close  examination  of  the 
CCT  value  data  (Table  21).  This  examination  revealed  that  four  sites, 

1,  3,  k,  and  2k,  exhibited  very  low  CCT  values  in  band  7.  An  examina- 
tion of  the  photomosaic  of  the  study  area  (Figure  33)  indicates  that 
sites  1,  3,  and  k are  all  in  floodplain  forest,  and  an  examination  of 
the  topographic  maps  (Figures  9 and  65)  reveals  that  those  three  sites 
were  probably  flooded.  On  the  other  hand,  site  2 is  also  in  a 
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floodplain  forest,  but  the  topographic  map  indicates  that  it  may  be  on  a 
low  ridge  so  that  it  was  probably  above  water  on  21  February  1975. 

136.  The  Mississippi  floodplain  forests  are  composed  almost 
exclusively  of  broadleaf  deciduous  species,  and  they  would  be  expected 
to  be  leafless  in  February.  Thus,  the  radiation  from  the  siirface  in 
floodplain  forest  areas  would  come  only  partly  from  the  trees;  some,  and 
perhaps  a major  component,  would  come  from  the  surface  beneath  the  trees. 
If  that  surface  were  water,  then  the  band-7  CCT  values  of  flooded 
forest  areas  would  be  expected  to  be  abnormally  low,  but  somewhat  higher 
th£in  open  water.  This  is  exactly  the  condition  suggested  by  the  CCT 
values  of  sites  1,  3,  and  2k  in  band  7. 

137.  This  led  to  the  possibility  that  open  water,  flooded  forest, 
and  dry  (i.e.  not  flooded)  forest  could  be  discriminated  by  a more  de- 
tailed subdivision  of  band-7  CCT  values.  With  this  in  mind,  the  cluster 
class  sets  in  Table  25  were  selected  for  testing.  Figure  66  shows  the 
resulting  map.  It  is  quite  apparent  that  many,  if  not  most,  of  the 
forested  areas  that  could  be  expected  to  be  flooded  are  correctly  clas- 
sified. The  scattering  of  small  patches  of  flooded  forest  mixed  among 
the  dry  forest  almost  all  occur  close  to  the  cirrved  border  between  the 
hilly  country  to  the  east  and  the  Mississippi  River  floodplain  to  the 
west.  Since  the  sxm  was  only  55  deg  above  the  horizon  to  the  southeast 
at  the  time  of  the  Landsat  overpass  (Table  2),  these  obviously  misclas- 
sified  areas  may  be  the  result  of  topographic  effects.  The  terrain 
tends  to  slope  steeply  to  the  northeast  along  the  valley  wall  (Figure  9), 
and  this  would  produce  very  low  CCT  values,  which  would  create  the  same 
general  effect  on  CCT  values  as  water,  especially  in  band  7. 

138.  An  examination  of  the  site  data  for  the  July  197^  overpass 
does  not  reveal  any  depression  of  band-7  CCT  values  for  sites  1,  3, 
and  h.  The  reason,  of  co\irse,  is  that  in  July  the  deciduous  trees  are 
in  full  leaf,  and  thus  virtually  all  of  the  energy  reaching  the  Landsat 
sensor  is  coming  from  the  forest  canopy  and  almost  none  from  the  forest 
floor . 

139»  However,  the  experience  with  band-7  data  in  the  February 
1975  data  suggested  that  anything  that  modified  the  open-water  surface 
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Study  area  showing  interpretation  of  water,  nonforest 
two  types  of  forest  from  21  February  1975  data 


might  be  detectable.  For  example,  by  July,  weeds  and  some  crops  are 
tall  enough  so  that  the  tops  might  project  above  the  water  surface,  pro- 
vided the  water  is  not  more  than  about  O.5  m deep.  If  so,  then  the  band- 
7 CCT  values  of  such  areas  woxild  be  somewhat  higher  than  those  repre- 
senting deep  water.  To  test  this  hypothesis,  the  cluster  class  limits 
of  cluster  class  set  1 (Table  20)  derived  from  the  July  197^  data  were 
altered  in  Table  26  by  dividing  the  water  range  in  band  7 into  two 
classes.  Figure  67  illustrates  the  consequences. 

lUO.  Since  the  forest  and  nonforest  cluster  class  limits  were  not 
altered,  these  two  classes  were  mapped  the  same  as  the  corresponding 
two  classes  in  Figirre  58.  However,  the  alteration  of  the  water  cluster 
class  limits  has  resulted  in  a separation  of  deep  and  shallow  water. 

The  deepwater  regions.  Dump  Lake  and  the  lowest  elevation  areas  (below 
the  85-ft  contour  lines  in  Figirre  9),  have  been  delineated  very  effec- 
tively, and  most  of  the  shallow-water  regions  (above  the  85-ft  contour) 
have  also  been  accurately  classified.  However,  the  Yazoo  River  is  clas- 
sified as  shallow  water,  which  it  obviously  is  not.  This  is  ijndoubtedly 
a misclassif ication  caused  by  spurious  signatures.  The  river  is  quite 
narrow  (Figure  25)  and  bordered  on  both  sides  by  bands  of  trees.  The 
river  is  not  quite  wide  enough  to  permit  a pixel  to  fit  into  the  water 
surface,  and  thus  all  of  the  pixels  falling  on  the  river  are  contaminated 
by  radiation  reflected  from  the  bordering  trees.  This  is  apparently 
enough  to  raise  the  band-7  CCT  values  above  the  limits  used  for  shallow 
and  deep  water  (Table  26). 
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Study  area  showing  interpretation  of  two  types  of  water 
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PART  IV:  CONCLUSIONS  AIR)  RECOMI-IENDATIONC 


Conclusions 


lUl.  Based  on  previous  experience,  it  may  be  concluded  that: 

a_.  The  interactive  procedure,  as  described  in  the  text  and 
outlined  in  Plates  1 and  2.  can  be  used  to  quickly 
obtain  criteria  by  which  certain  terrain  types  can  be 
reliably  classified  by  semiautomated  procedures. 

A set  of  criteria  derived  from  terrain  data  obtained  at 
one  season  cannot  be  used  to  classify  data  obtained  at 
another  season.  Both  spectral  reflects  \ce  characteristics 
and  absolute  colors  change  on  a seasonal  basis. 

c_.  Some  ground  truth  data  are  essential,  ince  they  are 
necessary  in  judging  the  acceptability  of  the  clusters 
of  CCT  values  used  as  the  basis  for  classification. 


Recommendations 

1U2.  Since  this  study  has  demonstrated  that  Landsat  spectral  data 
can  be  used  to  determine  areal  terrain  data  and  related  information  on 
in  situ  conditions  over  large  areas,  it  is  recommended  that  the  procedure 
be  further  demonstrated  and  assessed  in  a wide  variety  of  terrain  and 
environmental  conditions.  This  will  enhance  the  confidence  in  the  use 
of  the  overall  methodology. 

1U3.  Some  of  the  procedures  developed  under  this  study  have  not 
been  automated  to  permit  efficient  use.  As  a result,  it  is  recommended 
that  the  following  research  be  conducted  to  provide  an  automated  or 
semiautomated  method  of  interpreting  and  mapping  worldwide  terrain: 

a.  Develop  automated  procedures , and  write  a computer  pro- 
gram to  permit  more  efficient  extraction  of  radiance 
data  on  preselected  terrain  featvires  by  the  use  of 
multipass  Landsat  data  sets. 

Design  computer  procedures  to  analytically  compare  two 
sets  of  Landsat  data  for  the  same  region  for  different 
times  of  the  year  and  to  output  a map  portraying  the 
changes  in  the  conditions. 

£.  Develop  automated  procedures  for  selecting  gray-tone 

densities  fi^r  better  map  portrayal  of  terrain  and  land- 
use  conditions. 

lOU 
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Table  1 

Sub.lectively  Chosen  CCT  Value  Classes* 


CCT  Value 

CCT  Value  Ranges 

/Spectral  Band 

Class  No. 

Band  k 

Band  5 

Band  6 

1 

22-26 

15-18 

17-23 

5-11 

2 

26-28 

18-21 

23-26 

ll-lU 

3 

28-31 

21-23 

26-28 

ll*-l8 

k 

31-33 

23-26 

28-33 

18-23 

5 

33-35 

26-30 

33-35 

23-25 

6 

35-37 

30-33 

35-ltO 

25-31 

T 

37-‘+0 

33-36 

ii0-lt2 

31-39 

8 

36-Ui 

U2-U6 

39-^8 

9 

h6-^0 

I 

I 


* See  Figure  2. 


Characteristics  of  the  Three  Landsat  Spectral  Data  Sets  for  Area  Near  Vicksburg,  Mississippi 
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Table  3 

Calculated  Flight  Path  Angles  and  Pixel  Sizes 
for  Three  Landsat  Data  Sets 


®f 

°X 

Landsat 

Solution 

(Flight  Path  Angle) 

(Pixel  Width) 

(Pixel  Length) 

Data  Set 

Set* 

deg 

m 

m 

11  Jul 

1 

9.03 

58.3 

1974 

2 

9.13 

80.9 

Average 

9.08 

58.3 

80.9 

21  Feb 

1 

11.00 

57.5 

1975 

2 

9.13 

80.7 

Average 

10.07 

57.5 

80.7 

13  Oct 

1 

9.03 

58.3 

1975 

2 

8.59 

78.9 

Average 

8.81 

58.3 

78.9 

* Solution  sets  1 and  2 to  the  simultaneous  equations  exist  since  there 
were  three  unknowns  to  be  determined  from  four  equations. 


Table  U 


(Continued) 


• See  pareigraph  56  for  key  to  terrain  conditions. 


Table  U (Concluded) 


26 

27 

28 

29 

30 


11  Jul  1974 

A, 


Overflights 


21  Feb  1979 


13  0( 


"3 
A^ 


Table  5 


Signatures  from  Sample  Sites,  13  October  197^ 


Site 

Mean  and  Standard  Deviations,**  CCT 

Units 

No.* 

Band 

Band  5 

Band  6 

Band  7 

1 

15.8<0.7> 

14.8(1 .0) 

25. 6( 1 .0) 

14.0(0.9) 

2 

20.2(0.7) 

22.6(0.7) 

35.1(1.5) 

23. 1(1.5) 

3 

15.4(0.7) 

14.4(0.5) 

27. 4( 1 .9) 

16.6(1.7) 

4 

16.1(0.6) 

15.8(0.4) 

24.7(0.9) 

13.2(0.8) 

5 

21 .8(1 .2) 

25. 9( 1 .8) 

29.8(0.7) 

16.4(0. 7) 

6 

18.9(1.1) 

21 .6(0.9) 

32.2(2. 1 ) 

20.9(2.5) 

7 

20. 4( 1 .7) 

25.0(1.9) 

29.9(0.9) 

1 7.6( 1.5) 

8 

25.8(0.8) 

32.2(0.5) 

33. 0( 1 .2) 

16.9(0.6) 

9 

14.7(0.9) 

14.7(0.5) 

26. 8( 1.5) 

1 7.2(  1.6) 

10 

15.0(0.9) 

14.3(1 .0) 

23. 4( 1 .4) 

13. 9( 1.2) 

1 1 

13.9(0.6) 

14.0(0.9) 

24. 6( 1 .9) 

15.4(2.5) 

12 

1 5.6( 1 .0) 

14.3(0.7) 

24.9(0.8) 

14. 8( 1 .0) 

13 

16.3(0.5) 

16.0(1 .5) 

26.1(1.6) 

15.7(0.9) 

U 

17.6(2.3) 

20.7(3.9) 

28.3(2.4) 

16. 7(  1 . 7) 

15 

19.0(0.6) 

22. 2( 1.0) 

29.8(0.8) 

19.0(0.6) 

16 

17.8(1 .6) 

19.7(2.4) 

29. 2( 1.2) 

18.0(1.1) 

17 

21  .2( 1 .2) 

24.3(1 .3) 

32.0(2.2) 

19. 6( 1.6) 

18 

21.3(1.1) 

24.2(1.1) 

29.4(0.5) 

20. 7(0.5) 

19 

19.7(0.7) 

23.1(0.8) 

31.4(1.1) 

16.9(0.9) 

20 

23.4(0.9) 

28.4(0.7) 

30.9(0.6) 

16.0(0.7) 

21 

23.6(0.9) 

28. 7(0.9) 

31 . 7(0.7) 

1 7.3( 1.0) 

22 

21.1(1.1) 

26. 6( 1 .0) 

31.4(0.7) 

19.6(0.5) 

23 

25.6(0.5) 

33.6(0.7) 

35.2(0. 7) 

1 7. 8(0. 8) 

24 

18.1(1.2) 

16.6(1.2) 

12. 6( 1 .9) 

2.2( 1 .6) 

25 

25.2(0.4) 

32.1(0.7) 

34. 2( 1 .2) 

18.4(0.7) 

26 

19.7(0.7) 

24.6(2. 1 ) 

34. 2( 1 .0) 

22.3(0.9) 

27 

19.0(0.9) 

20.8(0.8) 

31 .8(0.7) 

21.3(0.7) 

28 

21.4(0.9) 

24.8(0.7) 

35.0(0.7) 

21 .9(0.8) 

29 

24.0(1 .3) 

28.6(2.6) 

32.1(1.5) 

1 7.  7(  1 . 7) 

30 

20.7(1 .0) 

24.6(0.5) 

31.8(1.1) 

19.3(0.5) 

* See 
**  Mean 

Figure  33  for  the 
and  standard  ( in 

locations  of 
parentheses) 

these  sites, 
deviations  are 

for  a 3-  by 

3-pixel  array. 


Sample 

Site  ±1.3  ±1.5  ±1.6  ±1.7  ±1.8 

26 


27  X X X X X 

28  X X X X X 

29 

30  _JL.  _JL  -JL  _JL  _JL_ 

Total  number 

meeting  criterion  15  17  19  19  20 


Table  7 

Effects  of  Variance  Window  on  CCT  Class  Limits 
13  October  1973  Data 


Class  Limits,*  CCT  Units 
Spectral  Band 


No.* 

Window 

Band  h 

Band  5 

Band  6 

Band  7 

1 

1.3 

15-17 

ll-l6 

21-27 

13-15 

1.5 

li*-17 

13-16 

21-27 

13-15 

1.8 

l4-l8 

13-17 

21-27 

12-16 

2 

1.3 

19-21 

21-21 

31-36 

22-21 

1.5 

19-22 

21-21 

31-37 

22-25 

1.8 

18-22 

21-21 

33-37 

21-25 

1* 

1.3 

15-17 

15-17 

23-26 

12-11 

1.5 

15-18 

11-17 

23-26 

12-15 

1.8 

l4-l8 

ll-l8 

23-26 

11-15 

5 

1.3 

21-23 

25-27 

29-31 

15-18 

1.5 

20-23 

21-27 

28-31 

15-18 

1.8 

20-21 

21-28 

28-32 

15-i8 

8 

1.3 

25-27 

31-33 

32-31 

l6-l8 

1.5 

21-27 

31-31 

32-31 

15-18 

1.8 

21-28 

30-31 

31-35 

15-19 

9 

1.3 

13-16 

13-16 

26-28 

l6-l8 

1.5 

13-i6 

13-16 

25-28 

i6-19 

1.8 

13-i6 

13-16 

25-29 

15-19 

10 

1.3 

ll-l6 

13-16 

22-25 

13-15 

1.5 

ll-l6 

13-16 

22-25 

12-15 

1.8 

13-17 

13-i6 

22-25 

12-16 

12 

1.3 

11-17 

13-i6 

21-26 

ll-l6 

1.5 

11-17 

13-16 

23-26 

13-16 

1.8 

11-17 

13-i6 

23-27 

13-17 

13 

1.3 

15-18 

15-17 

25-27 

11-17 

1.5 

15-18 

15-17 

25-28 

11-17 

1.8 

15-18 

ll-l8 

21-28 

ll-17 

15 

1.3 

18-20 

21-23 

29-31 

18-20 

1.5 

18-20 

21-21 

28-31 

18-20 

1.8 

17-21 

20-21 

28-32 

17-21 

18 

1.3 

20-23 

23-25 

28-31 

19-22 

1.5 

20-23 

23-26 

28-31 

19-22 

1.8 

20-23 

22-26 

28-31 

19-22 

(Continued) 

* Sites 

selected  using 

a screening 

window  of  ^1. 

8 CCT  units. 

Table  7 (Concluded) 


Site 

No. 

Variance 

Window 

Class  Limits,  CCT  Units 
Spectral  Band 

Band  U 

Band  5 

Band  6 

Band  7 

19 

1.3 

18-21 

22-24 

30-33 

16-18 

1.5 

18-21 

22-25 

30-33 

15-18 

1.8 

18-21 

21-25 

30-33 

15-19 

20 

1.3 

22-25 

27-30 

30-32 

15-17 

1.5 

22-25 

27-30 

29-32 

15-17 

1.8 

22-25 

27-30 

29-33 

l4-l8 

21 

1.3 

22-25 

27-30 

30-33 

16-19 

1.5 

22-25 

27-30 

30-33 

16-19 

1.8 

22-25 

27-30 

30-33 

16-19 

22 

1.3 

20-22 

25-28 

30-33 

18-21 

1.5 

20-23 

25-28 

30-33 

18-21 

1.8 

19-23 

25-28 

30-33 

18-21 

23 

1.3 

2I+-27 

32-35 

34-36 

17-19 

1.5 

2U-27 

32-35 

34-37 

16-19 

1.8 

2it-27 

32-35 

33-37 

16-20 

25 

1.3 

24-26 

31-33 

33-35 

17-20 

1.5 

24-27 

^1-34 

33-36 

17-20 

1.8 

23-27 

30-34 

32-36 

17-20 

27 

1.3 

18-20 

20-22 

31-33 

20-23 

1.5 

18-20 

19-22 

30-33 

20-23 

1.8 

17-21 

19-23 

30-34 

20-23 

28 

1.3 

20-23 

24-26 

34-36 

21-23 

1.5 

20-23 

23-26 

34-36 

20-23 

1.8 

20-23 

23-27 

33-37 

20-24 

30 

1.3 

19-22 

23-26 

31-33 

18-21 

1.5 

19-22 

23-26 

30-33 

18-21 

1.8 

19-22 

23-26 

30-34 

18-21 

Table  8 


Terrain  lYne  Assiainment  of  Sites  in  Cluster  Matrix* 


* See  Figure  h2. 


w 


Table  9 


Terrain  Types  Related  to  Site  Numbers, 
13  October  1975  Data 


Site  No. 


Terrain  Types 

Forest  Nonforest  Water 


Cluster 


1 

4 

10 

12 

2 

28 

8 

23 

25 

15 

18 

27 

30 

22 

19 

20 
21 


X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


I 

1 

I 

I 

II 

II 

III 
III 

III 

IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 


Used  in  Cluster  Matrix,  Figure  h2 
13  October  1975  Data 


Class  Limits,*  CCT  Units 
Spectral  Band 


No, 

Band  U 

Band  5 

Band  6 

Band  7 

Cluster 

1 

13-16 

21-27 

13-15 

1 

h 

15-18 

11-17 

23-26 

12-15 

10 

ll-l6 

13-16 

22-25 

12-15 

12 

Il-IT 

13-16 

23-26 

13-i6 

2 

19-22 

21-21 

31-37 

22-25 

II 

28 

20-23 

23-26 

31-36 

20-23 

8 

21-27 

31-31 

32-31 

15-18 

III 

23 

21-27 

32-35 

31-37 

16-19 

25 

21-27 

31-31 

33-36 

17-20 

15 

18-20 

21-21 

28-31 

18-20 

IV 

18 

20-23 

23-26 

28-31 

19-22 

27 

18-20 

19-22 

30-33 

20-23 

30 

19-22 

23-26 

30-33 

18-21 

22 

20-23 

25-28 

30-33 

18-21 

19 

18-21 

22-25 

30-33 

15-18 

20 

22-25 

27-30 

29-32 

15-17 

21 

22-25 

27-30 

30-33 

16-19 

I 


* Screening  window  = +1.5  CCT  units;  variance  window  = il.5  CCT  units. 


Table  11 

CCT  Classes  Defining  Spectral  Signatures  of  Sites 
Used  in  Cluster  Matrix,  Figure  Ul 
13  October  1973  Data 


Site 

Class  Limits,* 
Spectral 

CCT  Units 

Band 

No. 

Band  U 

Band  5 

Band  6 

Band  7 

1 

11+-17 

13-16 

24-27 

12-16 

2 

19-22 

21-24 

33-37 

21-25 

1* 

li*-l8 

14-17 

23-26 

12-15 

8 

2U-27 

31-34 

31-35 

15-19 

10 

13-17 

13-16 

22-25 

12-16 

12 

li+-17 

13-i6 

23-27 

13-16 

15 

17-21 

21-24 

28-31 

17-21 

18 

20-23 

23-26 

28-31 

19-22 

19 

18-21 

21-25 

30-33 

15-19 

20 

22-25 

27-30 

29-33 

14-18 

21 

22-25 

27-30 

30-33 

16-19 

22 

19-23 

25-28 

30-33 

18-21 

23 

2lt-27 

32-35 

34-37 

16-19 

25 

24-27 

30-34 

33-36 

17-20 

27 

17-21 

19-22 

30-33 

20-23 

28 

20-23 

23-26 

33-37 

20-24 

30 

19-22 

23-26 

30-33 

18-21 

* 


Screening  window  = +1.5  CCT  units;  variance  window  = ^.1.7  CCT  units. 


Table  12 


Clusters  and  Associated  CCT  Class  Ranges* 
13  October  1973  Data 


CCT  Value  Ranges 

Site  Spectral  Band 


Cluster** 

No. 

Band  U 

Band  5 

Band  6 

Band  7 

A 

2h 

16-20 

15-18 

ll-ll 

1-1 

Cluster  Class 

Limits 

16-20 

15-18 

ll-ll 

1-1 

B 

1 

i4-17 

13-16 

21-27 

12-16 

h 

11-18 

11-17 

23-26 

12-15 

10 

13-17 

13-16 

22-25 

12-16 

12 

lU-17 

13-16 

23-27 

13-16 

Cluster  Class 

Limits 

13-18 

13-17 

22-27 

12-16 

C 

2 

19-22 

21-21 

33-37 

21-25 

8 

21-27 

31-31 

31-35 

15-19 

15 

17-21 

21-21 

28-31 

17-21 

18 

20-23 

23-26 

28-31 

19-22 

19 

18-21 

21-25 

30-33 

15-19 

20 

22-25 

27-30 

29-33 

ll-l8 

21 

22-25 

27-30 

30-33 

16-19 

22 

19-23 

25-28 

30-33 

18-21 

23 

21-27 

32-35 

31-37 

16-19 

25 

21-27 

30-31 

33-36 

17-20 

27 

17-21 

19-22 

30-33 

20-23 

28 

20-23 

23-26 

33-37 

20-21 

30 

19-22 

23-26 

30-33 

18-21 

Cluster  Class 

Limits 

17-27 

19-35 

28-37 

1I-25 

* Variance  window  = +1.7  CCT  units. 

**  See  Figure  U4. 


Table  13 

Application  of  the  Cluster  Group  Stretching  Rule-.-  to  the 

Cluster  Separation  Logic  (A) | (b(c) 
T 3 


Cluster  Class  Limits,  CCT  Units 
Spectral  Band 


Cluster 

Band  U 

Band  5 

Band  6 

Band 

Before  Application 

of  Rules 

A 

16-20 

15-18 

11-1^4 

1-L 

B 

13-18 

13-17 

22-27 

12-16 

C 

1T-2T 

19-35 

28-37 

lL-25 

After  Omni group 

Rules 

A 

0-63 

15-18 

0-63 

I-I4 

B 

0-63 

13-17 

0-63 

12-16 

C 

0-63 

19-35 

0-63 

lL-25 

After  Intragroup 

Rules 

A 

0-63 

0-63 

0-63 

l-i+ 

B 

0-63 

0-18 

0-63 

12-16 

C 

0-63 

19-63 

0-63 

ll*-25 

After  Intergroup  Rules  1 and 

2 

A 

0-63 

0-63 

0-63 

0-i* 

B 

0-63 

0-18 

0-63 

12-16 

n 

0-63 

19-63 

0-63 

1^4-63 

After  Inter group 

Rule  3 

A 

0-63 

0-63 

0-63 

O-I4 

B 

0-63 

0-18 

0-63 

12-63 

C 

0-63 

19-63 

0-63 

12-63 

After  Intergroup 

Rule  1* 

A 

0-63 

0-63 

0-63 

0-8 

B 

0-63 

0-18 

0-63 

9-63 

C 

0-63 

19-63 

0-63 

9-63 

Table  lU 

Stretched  Cluster  Class  Sets,  13  October  1973 


Cluster  Class  Limits,  CCT  Units 
Spectral  Band 


Cluster 

Band  U 

Band  5 

Band  6 

Band  7 

Cluster  Class 

Set 

1 

A 

0-63 

0-63 

0-18 

0-63 

B 

0-63 

0-63 

19-27 

0-63 

C 

0-63 

0-63 

28-63 

0-63 

Cluster  Class 

Set 

2 

A 

0-63 

0-63 

0-63 

0-8 

B 

0-63 

0-18 

0-63 

9-63 

C 

0-63 

19-63 

0-63 

9-63 

Cluster  Class 

Set 

3 

A 

0-63 

0-63 

0-63 

0-8 

B 

0-63 

0-63 

0-27 

9-63 

C 

0-63 

0-63 

28-63 

9-63 

Cluster  Class 

Set 

1* 

A 

0-63 

0-63 

0-18 

0-63 

B 

0-63 

0-18 

19-63 

0-63 

C 

0-63 

19-63 

19-63 

0-63 

I 


Comparison  of  Areas  of  Classes,  Automated  Interpretations 
Versus  Air-Photo  Interpretation 


Table  l6 

Signatures  from  Sample  Sites,  11  July  197^ 


Site  Mean  and  Standard  Deviations,**  CCT  Units 


No. 

* Band  U 

Band  5 

Band  6 

Band  7 

1 

25.0(0.  > 

19.1(1.0) 

38.7(0.5) 

31 .4( 1.4) 

2 

25.4(0.9) 

19.0(0.9) 

38. 7(0.5) 

31  .8(1.3) 

3 

25.0(0.  > 

19.1(0.9) 

39.2(0.7) 

33. 3( 1 .0) 

4 

26. 0( 1.1) 

19. 2( 1 .0) 

40. 4( 1 .3) 

33.8(1.2) 

5 

32. 0( 1.7) 

30.6(1.3) 

41.0(0.9) 

29. 3( 1.3) 

6 

36.0(1 . 1 } 

35. 4( 1 .8) 

44.0(0.7) 

32. 3( 1.3) 

7 

30.8(1.5) 

29. 0( 1 .9) 

40. 7( 1 .0) 

28.4(1.8) 

8 

34. 8( 1 .5) 

35.1(1 .9) 

39. 2( 1 .3) 

24. 9( 1.1) 

9 

24.6(0.7) 

16.7(0.9) 

38.8(0.8) 

33.6(0.7) 

10 

25.0(0.  ) 

17.1(0.8) 

38.8(0.8) 

33. 0( 1 . 7) 

1 1 

24.7(0.5) 

16.3(0.5) 

37.2(0.8) 

31.1(2.2) 

12 

24.7(0.5) 

1 7.4( 1 .6) 

37. 8( 1 .6) 

31 .2(2.9) 

13 

24.9(0.9) 

1 7.2( 1.1) 

38.4(0.5) 

31.8(1.4) 

14 

27.2(2.5) 

22.9(4.1) 

37.7(0.9) 

29.1(1.3) 

15 

27. 8( 1 .7) 

24.3(2.5) 

36.2(1 .6) 

26. 6( 1 .9) 

16 

27. 2( 1.3) 

22.0(2.7) 

38. 6( 1 .5) 

30.1(1.6) 

1 7 

30.1(0.6) 

26.1(1.2) 

38.7(1 .0) 

28. 9( 1.6) 

18 

29.2(0.7) 

25. 0( 1.5) 

40.3(0.5) 

30. 7( 1.2) 

19 

27.7(0.9) 

22.0(0.9) 

39.3(1 .4) 

31.2(1.4) 

20 

28.6(0.9) 

24.4(0.7) 

41  .2(0.4) 

32. 2( 1.4) 

21 

30.8(0.7) 

27. 7( 1 .8) 

41 .9( 1.8) 

32.8(3.4) 

22 

26.8(0.7) 

20.1(0.6) 

42.3(0.9) 

35. 9( 1 .8) 

23 

31 .0(2.6) 

28.9(4.0) 

45. 8( 1 .0) 

39.2(3. 7) 

24 

26.3(1 .0) 

21 .2(0.7) 

20. 3( 1 .0) 

7. 7(0. 5) 

25 

26.7(0.9) 

19.9(0.6) 

48.4(0.7) 

49.9(2. 1 ) 

26 

31 .6(2.1 ) 

29.6(2.4) 

42. 4( 1 .6) 

32. 3( 1 .5) 

27 

31  .2(0.8) 

27. 2( 1 .4) 

47.7(0.5) 

41 .2(2.2) 

28 

29.5(0.7) 

27.9(4.5) 

43. 4( 1.1) 

34.7(3.5) 

29 

29. 3( 1 .4) 

26.2(2.1) 

42. 7(0.9) 

33.0(2.5) 

30 

31.3(2.3) 

30.1(3.2) 

43. 6( 1.1) 

33. 7( 1.4) 

See  Figure  33  Tor 

the  locations  of 

these  sites 

Mean  and  statidard 

(in  parentheses) 

deviations 

are  for  a 3- 

by  3-pixel  array. 
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Sample 

Site 


1 

2 

3 

1* 

5 

6 
T 
8 
9 

10 


11 

12 

13 

lU 

15 

16 

17 

18 

19 

20 


21 

22 

23 

2k 


25 

26 

27 

28 

29 

30 


1 


Table  17 

Effect  of  Changing  Screening  Window  Criterion  on  Number 


Total 

Accept- 

able 


of  Acceptable  Sample  Sites,  11  July  197^  and 
21  February  1975  Data 


X 


X 

X 


X 

X 


Screening  Window,  COT  Units 


+1.3 


+1.5 


+1.7 


X 

X 

X 


X 

X 


X 

X 

X 

X 


X 

X 


X 

X 

X 

X 

X 


X 

X 


X 

X 

X 


X 

X 


X 

X 

X 

X 


X 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 


X 

X 

X 

X 


X 

X 

X 

X 


X 

X 


X 

X 


10 


10 


12 


13  17 


+1.8 


+1.1 

IE  75  IE  H lE  H IE  H IE  75 


X 

X 

X 

X 

X 


X 

X 


X 

X 


X 

X 

X 

X 


15  17 


XX  XX  X X X X XX  XXXX 


Table  l8 

Effect  of  Variance  Window  on  CCT  Class  Limits 
11  July  197^  Data 


Site 

No.* 

Variance 

Window 

Class  Limits,  CCT  Units 
Spectral  Band 

Band  U 

Band  5 

Band  6 

Band  7 

1 

1.3 

2I+-26 

18-20 

37-40 

30-33 

1.5 

24-26 

18-21 

37-40 

30-33 

1.7 

23-27 

17-21 

37-40 

30-33 

1.8 

23-27 

17-21 

37-40 

30-33 

2 

1.3 

24-27 

18-20 

37-40 

31-33 

1.5 

24-27 

18-20 

37-40 

30-33 

l.T 

24-27 

17-21 

37-40 

30-33 

1.8 

24-27 

17-21 

37-40 

30-34 

3 

1.3 

24-26 

18-20 

38-40 

32-35 

1.5 

24-26 

18-21 

38-41 

32-35 

1.7 

23-27 

17-21 

38-41 

32-35 

1.8 

23-27 

17-21 

37-41 

32-35 

1.3 

25-27 

18-20 

39-42 

33-35 

1.5 

25-27 

18-21 

39-42 

32-35 

1.7 

24-28 

18-21 

39-42 

32-35 

1.8 

24-28 

17-21 

39-42 

32-36 

5 

1.3 

31-33 

29-32 

4o-42 

28-31 

1.5 

31-33 

29-32 

4o-42 

28-31 

1.7 

30-34 

29-32 

39-43 

28-31 

1.8 

30-34 

29-32 

39-43 

28-31 

6 

1.3 

35-37 

34-37 

43-45 

31-34 

1.5 

35-37 

34-37 

43-45 

31-34 

1.7 

34-38 

34-37 

42-46 

31-34 

1.8 

34-38 

34-37 

42-46 

31-34 

9 

1.3 

23-26 

15-18 

38-40 

32-35 

1.5 

23-26 

15-18 

37-40 

32-35 

1.7 

23-26 

15-18 

37-40 

32-35 

1.8 

23-26 

15-18 

37-41 

32-35 

10 

1.3 

24-26 

16-18 

38-40 

32-34 

1.5 

24-26 

16-19 

37-40 

32-34 

1.7 

23-27 

15-19 

37-40 

31-35 

1.8 

23-27 

15-19 

37-41 

31-35 

(Continued) 


Sites  selected  using  a screening  window  of  ^1.8  CCT  units. 


Table  l8  (Concluded) 


Class  Limits,  CCT  Units 


Site 

Variance 

Spectral 

Band 

No. 

Window 

Band  U 

Band  5 

Band  6 

Band  7 

13 

1.3 

2l*-26 

16-18 

37-i*0 

31-33 

1.5 

23-26 

16-19 

37-i*0 

30-33 

l.T 

23-27 

16-19 

37-**0 

30-33 

1.8 

23-27 

15-19 

37-i*0 

3O-3I+ 

IT 

1.3 

29-31 

25-27 

37-J+O 

28-30 

1.5 

29-32 

25-28 

37-1+0 

27-30 

1.7 

28-32 

2U-28 

37-1*0 

27-31 

X.8 

28-32 

2U-28 

37-1*0 

27-31 

18 

1.3 

38-30 

2lt-26 

39-1*2 

29-32 

1.5 

28-31 

2l*-26 

39-1*2 

29-32 

1.7 

28-31 

23-27 

39-1*2 

29-32 

1.8 

27-31 

23-27 

39-1*2 

29-32 

19 

1.3 

26-29 

21-23 

38-1+1 

30-32 

1.5 

26-29 

21-23 

38-1+1 

30-33 

1.7 

26-29 

20-2i* 

38-1+1 

30-33 

1.8 

26-29 

2O-2L 

38-1+1 

29-33 

20 

1.3 

27-30 

23-26 

1+0-1+2 

31-33 

1.5 

27-30 

23-26 

1+ 0-1+3 

31-3I+ 

1.7 

27-30 

23-26 

1+0-1+3 

31-31+ 

1.8 

27-30 

23-26 

39-1*3 

30-31+ 

22 

1.3 

26-28 

19-21 

1+1-1+ 1+ 

35-37 

1.5 

25-28 

19-22 

1+1-U1+ 

31+-37 

1.7 

25-28 

18-22 

1+1-1+ 1+ 

3I+-38 

1.8 

25-29 

18-22 

1+1-1+ 1+ 

3I+-38 

2h 

1.3 

25-28 

20-22 

19-22 

6-9 

1.5 

25-28 

20-23 

19-22 

6-9 

1.6 

25-28 

20-23 

19-22 

6-9 

1.8 

25-28 

19-23 

19-22 

6-9 

1 


Table  19 

CCT  Classes,  11  July  19TH  Data,  with  a 
Variance  Window  of  +1.3  CCT  Units 


Site 

Class 

Limits,  CCT  Units 
Spectral  Band 

No. 

Band  U 

Band  5 

Band  6 

Band  7 

2h 

25-28 

20-22 

19-22 

6-9 

1 

2J+-26 

18-20 

37-40 

30-33 

2 

2U-2T 

18-20 

37-40 

31-33 

3 

2U-26 

18-20 

38-40 

32-35 

h 

25-27 

18-20 

39-42 

33-35 

13 

24-26 

16-18 

37-40 

31-33 

9 

23-26 

15-18 

38-40 

32-35 

19 

26-29 

21-23 

38-41 

30-32 

20 

27-30 

23-26 

40-42 

31-33 

18 

28-30 

24-26 

39-42 

29-32 

Table  20 

Stretched 

Cluster  Class 

1 Sets,*  11  July  1974  Data 

CCT  Cluster  Class  Limits 
Spectral  Band 

Terrain 

Cluster 

Band  4 

Band  5 

Band  6 

Band  7 

Cluster 

Class  Set  1 

A 

0-63 

0-63 

0-63 

0-19 

Water 

B 

0-63 

0-20 

0-63 

20-63 

Forest 

C 

0-63 

21-63 

0-63 

20-63 

Nonforest 

Cluster 

Class  Set  2 

A 

0-63 

0-63 

0-29 

0-63 

Water 

B 

0-63 

0-20 

30-63 

0-63 

Forest 

C 

0-63 

21-63 

30-63 

0-63 

Nonforest 

As  eventually  used. 


Table  21 

SiKnatures  from  Sample  Sites,  21  P'ebruary  1973 


Site  Mean  and  Standard  Deviations,**  CCT  Units 


No.* ** 

Band 

Band  5 

Band  6 

Band  7 

1 

15.0(0.7) 

16.1(0.3) 

14.1(1.3) 

4.4(0. 5) 

2 

16.2(0.4) 

20.3(0.5) 

20.8(1.1) 

10.6(1.1) 

3 

14. 6( 1 .4) 

16.0(1.5) 

16.4(1.5) 

8.6(2. 5) 

4 

15. 3( 1 .2) 

16. 9( 1 .6) 

1 7.4( 1 .6) 

7.9( 1 .4) 

5 

16.3(0.5) 

17.1(1.1) 

38. 2( 1 .2) 

28.1(2.5) 

6 

16.6(0.5) 

19.9(0.6) 

28.4(0.5) 

16. 0( 1.2) 

7 

17.3(1 .2) 

21 .7( 1 .4) 

23. 8( 1 .0) 

1 1 .0(0.9) 

8 

16.9(0.6) 

19.6(0.5) 

25.1(0.8) 

12.7(0.9) 

9 

13.4(1.1) 

15.3(0.9) 

19. 6(  1 . 7) 

10. 9( 1 . 5) 

10 

13.1(0.8) 

14. 6( 1.1) 

21.6(1 .6) 

10.6(2.2) 

1 1 

13.3(0.5) 

15. 8( 1.2) 

21 .3( 1 .0) 

1 1 . 7(0.9) 

12 

13.3(0.9) 

14. 3( 1 .0) 

22.6(3.2) 

12.0(2.0) 

13 

14.9(0.8) 

1 7. 7( 1 .3) 

23. 7( 1 .6) 

13.0(0. 7) 

14 

15. 6( 1.1) 

17. 0( 1 .9) 

21 .7(2.9) 

14.9(4.0) 

15 

16.3(2.5) 

18.1(3.0) 

28.2(4.3) 

21  .2(1.3) 

16 

17. 2( 1.1) 

21  .4( 1 .7) 

30. 8(  1.7) 

21.0(2.4) 

17 

17.1 (0.6) 

19.1(0.6) 

31 .1(0.9) 

19.3(0. 7) 

18 

17.2(0.4) 

20.0(0.  ) 

28.2(0.8) 

15.3(0.5) 

19 

17.4(0.7) 

20.8(0.8) 

25. 8( 1 .2) 

1 3.2( 1.0) 

20 

17.7(0.7) 

21 .3( 1 .0) 

28.9(0.8) 

1 7.3( 1 . 1 ) 

21 

22. 7( 1 .0) 

28. 8( 1 .6) 

32.4(3.4) 

20.6(2.  1 ) 

22 

17.0(0.7) 

21 .6(0.7) 

26.6(0.7) 

13.8(1.1) 

23 

19.5(1 .0) 

22. 4( 1 . 7) 

33. 8( 1 .4) 

20.0(0.9) 

24 

20. 3( 1.8) 

24.4(0.5) 

16.9(1 .2) 

2. 2(2.0) 

25 

17.4(0.7) 

20.1 (0.6) 

32. 7(3.4) 

22. 9( 1 .5) 

26 

21 .0(0.9) 

24. 9( 1 .3) 

32.4(2.5) 

18. 2( 1 .0) 

27 

22.2(0.7) 

26. 3( 1.7) 

29.9(2.2) 

1 7.3(0. 7) 

28 

18.2(0.7) 

22.0(0.9) 

28.8(0.8) 

16. 6( 1.2) 

29 

22. 5( 1.2) 

28. 8( 1 .6) 

31  .0(1.3) 

18.2(1.1) 

30 

21.4(0.7) 

24.2(2.2) 

32.1(2.6) 

20. 2( 1.2) 

* See  Figure  33  for  the  locations  of  these  sites, 

**  Mean  and  standard  (in  parentheses)  deviations  are  for  a 3- 
by  3-pixel  array. 
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Table  22 


Effect  of 

Variance  Window 

on  CCT  Class 

Limits 

21  February  1975  Data 

Site 

Variance 

Class  Limits,*  CCT  Units 
Spectral  Band 

No.* 

Window 

Band  U 

Band  5 

Band  6 

Band  7 

1 

1.3 

lit-l6 

15-17 

13-15 

3-6 

1.5 

lh-l6 

15-18 

13-16 

3-6 

1.8 

13-17 

lit-18 

12-16 

3-6 

2 

1.3 

15-17 

19-22 

20-22 

9-12 

1.5 

15-18 

19-22 

19-22 

9-12 

1.8 

lU-l8 

19-22 

19-23 

9-12 

k 

1.3 

i4-17 

16-18 

16-19 

7-9 

1.5 

i4-17 

15-18 

16-19 

6-9 

1.8 

lU-17 

15-19 

16-19 

6-10 

6 

1.3 

15-18 

19-21 

27-30 

15-17 

1.5 

15-18 

18-21 

27-30 

15-17 

1.8 

15-18 

18-22 

27-30 

1I+-18 

7 

1.3 

16-19 

20-23 

23-25 

10-12 

1.5 

16-19 

20-23 

22-25 

10-12 

1.8 

16-19 

20-23 

22-26 

9-13 

8 

1.3 

16-18 

18-21 

2I-26 

II-II+ 

1.5 

15-18 

18-21 

2U-27 

11-li* 

1.8 

15-19 

18-21 

23-27 

11-14 

9 

1.3 

12-15 

1I-I7 

18-21 

10-12 

1.5 

12-15 

lU-17 

18-21 

9-12 

1.8 

12-15 

lit-17 

18-21 

9-13 

11 

1.3 

12-15 

15-17 

20-23 

10-13 

1.5 

12-15 

1I-I7 

20-23 

10-13 

1.8 

12-15 

li+-l8 

20-23 

10-13 

13 

1.3 

lU-16 

16-19 

22-25 

12-14 

1.5 

13-16 

16-19 

22-25 

12-14 

1.8 

13-17 

16-19 

22-25 

11-15 

IT 

1.3 

16-18 

18-20 

30-32 

18-21 

1.5 

16-19 

18-21 

30-33 

18-21 

1.8 

15-19 

17-21 

29-33 

18-21 

(Continued) 


* Sites  selected  using  a screening  window  of  j^l.8  CCT  units. 
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Table  22  ( Concluded ) 


Site 

Variance 

Class  Limits,  CCT  Units 
Spectral  Band 

No. 

Window 

Band  h 

Band  5 

Band  6 

Band  7 

18 

1.3 

16-18 

19-21 

27-29 

14-17 

1.5 

16-19 

19-21 

27-30 

14-17 

1.8 

15-19 

18-22 

26-30 

14-17 

19 

1.3 

16-19 

20-22 

25-27 

12-14 

1.5 

16-19 

19-22 

24-27 

12-15 

1.8 

16-19 

19-23 

24-28 

11-15 

20 

1.3 

16-19 

20-23 

28-30 

16-19 

1.5 

16-19 

20-23 

27-30 

16-19 

1.8 

16-19 

20-23 

27-31 

16-19 

22 

1.3 

16-18 

20-23 

25-28 

13-15 

1.5 

16-18 

20-23 

25-28 

12-15 

1.8 

15-19 

20-23 

25-28 

12-16 

23 

1.3 

18-21 

2I-2U 

33-35 

19-21 

1.5 

18-21 

21-24 

32-35 

19-21 

1.8 

18-21 

21-24 

32-36 

18-22 

28 

1.3 

iT-19 

21-23 

28-30 

15-18 

1.5 

17-20 

21-23 

27-30 

15-18 

1.8 

16-20 

20-24 

27-31 

15-18 

29 

1.3 

21-2I* 

28-30 

30-32 

17-19 

1.5 

21-2U 

27-30 

30-32 

17-20 

1.8 

21-2k 

27-31 

29-33 

16-20 

r 


Table  23 

CCT  Classes.  21  February  1975  Data,  with  a 
Variance  Window  of  +1.5  CCT  Units 


Class  Limits,  CCT  Units 


Site  Spectral  Band 


No. 

Band  k 

Band  5 

Band  6 

Band 

2h 

19-22 

23-26 

15-18 

1-4 

1 

1I+-I6 

15-18 

13-16 

3-6 

11 

12-15 

li+-17 

20-23 

10-13 

2 

15-18 

19-22 

19-22 

9-12 

7 

16-19 

20-23 

22-25 

10-12 

8 

15-18 

18-21 

2l*-27 

ll-l4 

19 

16-19 

19-22 

2U-2T 

12-15 

22 

16-18 

20-23 

25-28 

12-15 

18 

16-19 

19-21 

27-30 

i4-17 

6 

15-18 

18-21 

27-30 

15-17 

28 

17-20 

21-23 

27-30 

15-18 

20 

16-19 

20-23 

27-30 

16-19 

IT 

16-19 

18-21 

30-33 

18-21 

Table  2U 

Stretched  Cluster  Class  Sets,*  21  February  1975  Data 


Cluster  Class  Limits,  CCT  Units 
Spectral  Band 

Cluster 

Band  4 

Band  5 

Cluster  Class  Set  1 

Band  6 

Band  7 

A 

0-63 

21-63 

0-18 

0-63 

®2 

0-63 

0-20 

0-18 

0-63 

0-63 

0-17 

19-63 

0-63 

C 

0-63 

18-63 

Cluster  Class  Set  2 

19-63 

0-63 

A 

0-63 

21-63 

0-63 

0-7 

B, 

0-63 

0-20 

0-63 

0-7 

B^ 

2 

0-63 

0-17 

0-63 

8-63 

C 

0-63 

18-63 

0-63 

8-63 

* With  two  types  of  forests,  and 


I 


J 


r 


Table  25 

Cluster  Class  Sets  Selected  to  Discriminate  Between 
Flooded  and  Unflooded  Forest 
21  February  1975  Data 


Cluster  Class  Limits,  CCT 
Spectral  Band 

Units 

Terrain 

Cluster 

Band  ^4 

Band  5 

Band  6 

Band  7 

Type 

A 

0-63 

21-63 

0-63 

0-7 

Open  water 

0-63 

0-20 

0-63 

0-7 

Flooded  forest 

®2 

0-63 

0-17 

0-63 

8-63 

Dry  forest 

C 

0-63 

18-63 

0-63 

8-63 

Nonforest 

Table  26 

Revised  Stretched  Cluster  Class  Set  1,  11  July  ig'jh  Data 


Cluster  Class  Limits,  CCT 
Spectral  Band 

Units 

Terrain 

Cluster 

Band  L 

Band  5 

Band  6 

Band  7 

Type 

0-63 

0-63 

0-63 

0-8 

Deep  water 

^2 

0-63 

0-63 

0-63 

9-19 

Shallow  water 

B 

0-63 

0-20 

0-63 

20-63 

Forest 

C 

0-63 

21-63 

0-63 

20-63 

Nonforest 

PROCEDURE  FOR  SELECTING 
APPROPRIATE  CLASSIFICATION  CRITERIA 


APPENDIX  A:  DESCRIPTION  OF  OPTRONICS  FILM  READER/WRITER 


1.  Maps  depicting  the  interpi-eted  terrain  conditions  ai'e  produced 
using  an  incremental  film  reader/writer  (Figure  Al).  The  film  reader/ 
wi'iter  is  an  electromechanical  photograph-scanning  and  film-writing 


Figure  Al.  Film  reader/writer 


system  designed  to  produce  images  up  to  22.3  by  cm  in  size.  The 

instrument  can  operate  in  either  of  two  modes — an  input  or  scanning  mode, 
or  an  output  or  writing  mode.  Tlie  writing  mode  is  useti  to  convert  CCT 
values  to  images  in  pliotographic  foi'in.  Foi’  tliis  mode  of  operation,  tiie 
instrument  i.3  equipped  wLtii  a rotating  drum  and  an  optical  .system  con- 
sisting of  a light-emitting  diode  (LED),  a selectable  ajierture,  and  a 
lens  system  that  focuses  a spot  of  light  from  the  LED  onto  tiie  perimeter 
of  tiie  drum.  Tiie  drum  is  lioused  in  a 1 ight-t  i g.iit  enclosure,  wiiicli  is 
demountaliLe  and  is  removed  to  a piiotograplii c dai'kj’oom  for  film  ioailinr 
and  unloading.  A piece  of  film  is  cirunped  to  tiie  outside  of  tiie  drum 
for  exposure. 

2.  Til'.'  film  can  in;  exposed  us.ing  12.h-i  or  'jO.O-iim  siiots 


Ai 


of  light  from  the  LED.  The  raster  interval  and  spot  size,  which  are 
selectable,  control  the  size  of  the  negatives  that  will  be  produced.  As 
the  drum  rotates,  the  carriage  supporting  the  optical  system  is  stepped 
one  step  per  drum  revolution  in  the  axial  direction  at  the  selected 
raster  interval  until  the  total  area  of  the  film  or  the  area  of  interest 
has  been  exposed.  The  use  of  high-speed  film  permits  very  short  ex- 
posure times  and  results  in  a recording  rate  of  up  to  60,000  exposures 
(or  pixels)  per  second. 

3.  The  intensity  of  light  from  the  LED  is  modulated  incrementally 
in  proportion  to  the  values  recorded  on  the  magnetic  tape.  Thus,  as  the 
drum  of  the  film  reader /writer  rotates,  a spot  is  exposed  on  the  film 
for  each  pixel  value  in  a scan  line  (row).  When  an  interrecord  gap  code 
occurs,  the  LED  is  extinguished  until  the  drum  revolution  is  completed 
and  the  carriage  for  the  optical  system  has  moved  forward  one  increment. 
Exposure  of  the  next  row  of  spots  then  begins.  This  process  is  repeated 
until  the  end-of-file  code  occurs  on  the  magnetic  tape. 

h.  The  film  reader/writer  is  controlled  by  a minicomputer  so  that 
real-time  manip\ilation  of  the  CCT  values  can  be  used  to  produce  a number 
of  photographic  effects.  However,  for  this  capability  to  be  used 
successfully,  very  careful  control  over  the  photographic  process  and 
knowledge  of  the  relation  between  the  digital  input  from  the  magnetic 
tape  and  the  photographic  output  of  the  film  writer  are  required. 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimilfe  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Struve,  Horton 

Acquisition  of  terrain  information  using  Landsat  multispectral 
data;  Report  2:  An  interactive  procedure  for  classifying  terrain 
types  by  spectral  characteristics  / by  Horton  Struve,  Warren  E. 
Grabau,  and  Harold  W.  West.  Vicksburg,  Miss.  : U.  S.  Waterways 
Experiment  Station,  1977. 

105,  [28j,  2 p.,  2 leaves  of  plates  : ill.  ; 27  cm.  (Tech- 
nical report  - U.  S.  Army  Engineer  Waterways  Experiment  Station  ; 
M-77-2,  Report  2) 

Prepared  for  Assistant  Secretary  of  the  Army  (R&D) , Department 
of  the  Army,  Washington,  D.  C.,  under  Project  4A161101A91D, 

Task  02,  Work  Unit  095. 
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